REPORT 976 


LINEAR THEORY OF BOUNDARY EFFECTS IN OPEN WIND TUNNELS WITH 

FINITE JET LENGTHS 


Bv S. Katzoff, Clifford S. Gardner, Leo Diesendruck, and Bertram J. Eisenstadt 


SUMMARY 

In the first part, the boundary conditions for an open wind 
tunnel (incompressible flow) are examined with special refer- 
ence to the effects of the closed entrance and exit sections. Basic 
conditions are that the velocity must be continuous at the en- 
trance lip and that the velocities in the upstream and down- 
stream closed portions must be equal. For the two-dimensional 
open tunnel, interesting possibilities develop from the fact that 
the pressures on the two free surfaces need not be equal. 

Electrical analogies that might be used for solving the flow in 
open wind tunnels are outlined. Two types are described — 
one in which electrical potential corresponds to velocity poten- 
tial, and another in which electrical potential corresponds to 
acceleration potential. The acceleration-potential analogies 
are probably experimentally simpler than the velocity-potential 
analogies. 

In the second part, solutions are derived for four types of 
two-dimensional open tunnels, including one in which the 
pressures on the two free surfaces are not equal. Numerical 
results are given for every case. In general, if the lifting ele- 
ment is more than half the tunnel height from the inlet, the 
boundary effect at the lifting element is the same as for an 
infinitely long open tunnel. 

In the third part, a general method is given for calculating 
the boundary effect in an open circular wind tunnel of finite 
iet length. Numerical results are given tor a lifting element 
concentrated at a point on the axis. 

INTRODUCTION 

The basic theory of boundary corrections for an open 
wind tunnel was given by Prandtl many years ago (reference 
1) and has since been used with reasonable success. The 
infinitely long open jet that was assumed in Prandtl’s 
analysis, however, has been frequently questioned as an 
adequate representation for an open wind tunnel, which 
normally has a relatively short jet between closed entrance 
and exit regions. The present examination of the problem 
was occasioned by the need for boundary corrections for 
tests in the Langley full-scale tunnel of a large helicopter, of 
which the forward edge of the rotor disk reached almost to 
the mouth of the entrance bell while the rear edge approached 
the exit bell. Previous studies (reference 2) had shown that 
the Prandtl theory was satisfactory for a wing in the usual 
position in the tunnel (about 20 feet downstream of the 


entrance), but it was felt that this simple theory was in 
adequate for such far forward and rearward locations of the 
lifting surface and that some further development was 
desirable. The only previous analysis bearing directly on 
the problem seemed to be that of reference 3, which con 
sidered a lifting element concentrated at a point on the axis 
of a circular open tunnel of finite jet length; however, the 
treatment therein was not rigorous and was justified only by 
a somewhat heuristic discussion, so that its general applica- 
bility was not obvious. Other studies treated either two- 
dimensional or axially symmetrical conditions (references 4 
and 5) and also did not consider the closed exit region, so 
that the extent of their applicability to the present problem 
was not at first apparent. A similar German wartime report 
(reference 6), which did not become available until after the 
present report was written, would have been more useful in 
this respect because of the generality of its physical discussion. 

Because of the particular shape of the tunnel cross section 
a reasonably simple solution in terms of available functions 
seemed unlikely; accordingly, the initial effort was directed 
toward defining the problem in such a way that it could be 
solved by analogy methods in an electrical tank. Identifica- 
tion of the necessary boundary conditions appeared at first 
to be somewhat perplexing; however, after recognition of 
some of the basic physical phenomena, the boundary condi- 
tions were readily clarified. The problem is thus considered 
now to be fairly well understood, at least insofar as it can be 
considered linear and uninfluenced by turbulent mixing at 
the free surfaces or by the irregular nature of the flow at the 
exit. As will appear later, however, grave technical difficul- 
ties exist in the exact solution by electrical-analogy methods 
so that, for example, rigorous evaluation of the tunnel inter- 
ference for the large helicopter in the Langley full-scale 
t unn el, which problem instigated the present research, has 
not yet been attempted. 

After the boundary conditions were clarified, analytical 
methods of solution were developed for two-dimensional and 
circular open tunnels. These studies have been combined 
with the discussion of the boundary conditions and the 
electrical analogies to form the present report, which, it is 
hoped, will serve to clarify basic concepts and establish a 
sound basis for any further work. • 

The report is divided into three parts. In part I, the 
boundary conditions are defined and discussed for the open 
wind tunnel with closed entrance and exit sections, and an 
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outline is given of suggested electrical analogies applicable 
to tbe problem. In part II, analytical solutions are given 
for various two-dimensional open-tuniiel types, together 
with numerous calculated results. In part III, a method of 
solution for the circular open tunnel is given, together with 
numerical results for the case of a lifting element concen- 
trated at a point on the axis of the tunnel. The treatment 
in every case is a linear one in which deformation of the jet 
boundary is considered to be small. , 

The parts were essentially independently prepared. 
Messrs. Gardner and Diesendruck contributed the analysis 
of part II. Mr. Eisenstadt contributed part III. Dr. 
Katzoff contributed part I and, in the absence of the others, 
prepared the numerical results of part II, made several 
minor revisions, and served as general editor of the report. 

I. BOUNDARY CONDITIONS AND ELECTRICAL 
ANALOGIES 

In part I, boundary conditions for an open wind tunnel 
are discussed with special reference to the effects of the 
closed entrance and exit sections. It is shown that the 
velocity on the free surface is not necessarily equal to the 
velocity far upstream in the closed portion and that cross- 
flows may exist in the free surface, unlike the case of the 
infinitely long open jet. A basic condition — analogous to 
the Kutta-Joukowski condition for the flow at the trailing 
edge of an airfoil — is that the velocity be continuous at the 
entrance lip. Electrical analogies that might be used for 
solving the flow in open wind tunnels are outlined. Two 
types are described; — one in which electrical potential corre- 
sponds to velocity potential, and another in which electrical 
potential corresponds to acceleration potential. 

BOUNDARY CONDITIONS 

Resume of Prandtl’s theory. — In PrandtI’s original dis- 
cussion, in which the entrance and exit regions are neglected, 
the open jet is considered as an infinitely long cylinder on the 
entire surface of which the pressure is constant, whence, 
by Bernoulli’s law, the velocity on the surface is constant. 
If this velocity is considered as the sum of the undisturbed 
tunnel velocity V and a small perturbation velocity (u,v,w) 
due to the presence of a body in the jet, the condition is then that 

(27 + ii) 2 + « 2 +w 2 «C/ 2 +227ii=Constant 

from which it is concluded that u is constant over the 
entire surface. Furthermore, since u is obviously zero far in 
front of the body, it must be zero over the entire surface. 

A corollary is that, on the jet surface, the perturbation 
velocity normal to u (that is, the circumferential velocity) 
is also zero, as is readily shown from a consideration of the 
rectangular path SPQR on the surface of the jet. (See fig. 
1 (a).) As has just been shown, the velocity component u 
parallel to the lines SP and QR is zero; hence, the perturba- 
tion potentials at points P and Q are the same as at points S 




Figure 1.— Illustrations for discussion of surface perturbation velocity in open wind tunnels. 

and R, respectively. If points S and R are far upstream 
of the body their potentials will be equal, so that the po- 
tentials at points P and Q are equal. The perturbation 
potential is thus uniform over the entire surface, and only 
perturbation velocities normal to the surface can exist at 
the surface. 

Modification of basic concepts. — If the closed entrance 
region is near the body, as shown in figure 1 (b), the preceding 
discussion and conclusions no longer apply. Thus, although 
u must still be constant ove^ the entire free surface, it is no 
longer necessarily zero; that is, the total velocity on the free 
surface is not necessarily equal to the velocity far upstream 
in the closed portion of the tunnel. The two velocities will, 
in fact, generally be unequal except in special cases where 
equality results from geometrical symmetry of the arrange- 
ment. (For example, if a horseshoe vortex is located in the 
horizontal plane of symmetry of the tunnel, the values of u 
at the top and bottom of the tunnel would be expected to be 
equal and opposite; but since u must be uniform over the 
surface, it follows that «=0.) Furthermore, the velocities 
in the jet surface normal to u (that is, the circumferential 
velocities) are, in general, no longer zero (except for axially 
symmetrical flows, such as that produced by a body of 
revolution on the axis of a circular tunnel) so that two 
surface points at the same longitudinal position, as P' and 
Q' (fig. 1 (b)), do not necessarily have the same values of tho 
perturbation potential. 
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Entrance-lip condition. — Consider, for simplicity, the 
sy mm etrical case of figure 1 (b), in which the lifting element 
is on the horizontal plane of symmetry of the tunnel. Since 
u is zero on the free boundary, the perturbation potential is 
constant, along the elements AB, CD, EF, . . ., although, 
as just indicated, it is not necessarily the same for all these 
elements. This one boundary condition for the open sec- 
tion — that the potential be constant along each of these 
elements — does not suffice, however, to define the problem 
uniquely. In fact, as will be obvious from the subsequent 
discussion of electrical analogies, the potentials of these 
elements may be quite arbitrarily assigned without violating 
this condition or the boundary condition on the closed por- 
tion of the tunnel (that the normal derivative of the potential 
be zero at the wall) . In order to avoid this lack of unique- 
ness, further conditions must be sought . The most important 
of these is that the velocity be continuous (in particular, 
not infinite) at the entrance lip (points A, C, E, . . .). 
This condition takes cognizance of the fact that, because of 
viscosity, the physical flow leaves the lip smoothly, just as 
it leaves the trailing edge of an airfoil; the condition is, in 
fact, strictly analogous to the Kutta-Joukowski condition 
for the trailing-edge of an airfoil, which similarly takes into 
account the basic viscosity effect and provides uniqueness 
where otherwise an infinity of solutions would exist. It is 
recognized that, just as the Kutta-Joukowski condition does 
not always suffice to predict airfoil lift very accurately, the 
corresponding condition for the open tunnel may similarly 
oversimplify the entrance-lip flow; however, as with the 
airfoil, the condition is probably adequate where the flow 
is not subject to an excessive pressure rise on approaching 
the lip. References 4, 5, and 6 used the condition, and 
reference 6, in addition, discussed it from the physical view- 
point and compared it with the airfoil trailing-edge condition. 

Concerning the downstream end of the open section, the 
exit lip may be considered to correspond to the leading edge 
of an airfoil and no effort need be made in an idealized flow 
analysis to eliminate infinite values of u at this edge. 

Jet contraction or expansion. — It has already been pointed 
out that, with a body in the jet, the velocity on the free 
surface is not necessarily equal to the velocity far upstream 
in the closed portion. During the course of the investiga- 
tion, it was noted that solutions could be obtained showing 
a difference between these two velocities, even when there 
was no body in the jet. Such a flow corresponds merely to 
a contraction or expansion of the jet, as indicated in figure 2. 
Thus, in figure 2 (a) , the velocity on the free surface is lower 
than the upstream velocity and remains so even as it ap- 
proaches the exit, in spite of the gradual contraction of the 
jet, because of the continuously increasing surface curvature. 
The velocity suddenly increases at the exit lip and finally 
is established at a value greater than that of the upstream 
velocity. With reasonable ratios of entrance to exit area, 
the flows of figure 2 may be readily obtained experimentally. 



(b) 

(a) Contracting jet. The pressure on the free surfoce exceeds both, the upstream and down- 

stream pressures but is very close to the upstream pressure. 

(b) Expanding Jet. The pressure on the free surface is less than either the upstream or down- 

stream pressure but is very close to the upstream pressure. 

Figure 2.— Contracting and expanding jets (2 or 3 dimensions). 

The significance of this expanding or contracting flow is 
that it represents a solution that satisfies all the boundary 
conditions previously discussed and is nevertheless undesir- 
able. In order to avoid such solutions, a further condition 
must accordingly be recognized; namely, that the velocities 
in the closed portions far upstream and far downstream of 
the open section be equal. 

It may be objected that in the normal design of an open 
wind tunnel the exit section is made larger than the entrance 
section. The purpose of the increased area is to allow for 
the reduced velocity toward the surface of the jet resulting 
from turbulent mixing with the surrounding still air. In- 
creasing the exit area by other than the correct amount will 
result in the type of flow indicated in figures 2 (a) or 2 (b) 
with a corresponding velocity gradient along the center of 
the t unn el. In any potential-flow solution these viscous 
effects cannot be considered. 

Spillage. — When an airfoil is tested at a high lift coefficient 
in an open tunnel, the downward deflection of the jet may 
result in appreciable spillage from the lower lip of the exit, 
together with lack of contact of the main flow with the upper 
lip. (See fig. 3 (a).) The air lost by spillage is replaced by 
air (of, however, a lower total pressure) entrained in the 
exit. Even without otherwise considering the distortion of 
the free surface, these Sow characteristics might seem too 
much at variance with the previously assumed characteristics 
to permit application of the theories being discussed. The 
calculations of part II for the two-dimensional open tunnel 
(that is, a rectangular tunnel with closed sides but open top 
and bottom) show, however, very little difference between the 
tunnel-induced-downwash distributions for the tunnel with 
two exit lips and the tunnel with one exit lip. That is, if 
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figure 3 (a) is assumed to represent a two-dimensional. flow, 
the fact that the upper lip of the exit is out of the flow field 
so that the lower lip takes over the entire burden of straight- 
ening the jet does not greatly affect the induced downwash. 

The effect of the exit lip on the flow phenomena is the least 
clear of the various phases of the present problem. For open 
wind t unnel s having essentially unflared exits, similar to that 
indicated in figure 3 (a), the suggestions of the preceding 
paragraph are probably adequate. The exit of the Langley 
full-scale tunnel, however, has a large bell mouth, and when 



(a) 






(b) 




(o) 


(a) Straight exit. 

(b) Bell-mouth exit. 

(c) Enclosed space beneath the lower tree surface (two-dimensional tunnel). 
Figure 3.— Spillage from the lower lip of the exit. 


airplanes are being tested at high lift coefficients a downward 
deflection of the air off the lower part of the bell, roughly 
as indicated in figure 3 (b), occurs. Whether the previously 
suggested concepts or, indeed, any linear theory can serve 
satisfactorily for this case seems questionable. 

Unequal surface pressures. — Aninterestingmethod of avoid- 
ing spillage suggests itself in the case of the two-dimensional 
open tunnel: If the space below the tunnel is enclosed, an 
excess pressure will be built up in this space, compared 
■with the pressure in the space above the upper free surface, 
so that the flow will be pushed up sufficiently to eliminate 
the spillage and insure precise contact of the lower free 
surface with the lower exit lip. (See fig. 3 (c).) The extent 
to which a free two-dimensional jet can be deformed bj r a 
pressure difference across its boundaries or, stated differently, 
the extent to which a two-dimensional free jet will deform in 
order to follow the only available path, is indicated by the 
smoke-flow photograph in figure 4. The setup consisted 
merely of a two- d imensional open jet -with entrance and exit 
sections displaced vertically relative to each other, arranged 
between transparent side walls, and provided with enclosed 
spaces above and below. 



Figure 4.— Two-dimensional jet with different pressures on the two free surfaces. 
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Details of interest in the figure, in addition to the jet 
deformation, are: 

(1) Separation of the flow from the upper lip of the exit 
because of the large angle of entry. A small bell mouth at 
the exit lip might have prevented such separation. 

(2) The rough flow on the upper surface at the entrance, 
compared with the smooth flow on the lower surface, reflects 
the fact that the boundary layer approaching the entrance 
is subjected to a rising pressure on the upper surface and a 
dropping pressure on the lower surface. 

(3) Because of turbulent mixing at the free surfaces, a 
certain amount of the air in the closed chambers above and 
below the jet is entrained in the jet. An equivalent quantity 
must be released, or skimmed off, at the exit in order that the 
total quantity in each chamber remain constant. This cir- 
culating mechanism results in the apparent overflow at the 
two exit lips. The return of the s kimm ed -off part to the jet 
surface can be seen at the bottom of the photograph. 

Tunnel without a closed exit. — Some mention is made in 
the subsequent discussion of the hypothetical open tunnel 
having a closed upstream entrance region but no closed exit 
region, the open section thus extending downstream to 
infinity. Calculations for such an arrangement (see part II) 
are generally simpler than for the actual tunnel with the 
closed exit and give very nearly the same answer, provided 
that the region of interest is much closer to the entrance 
than to the exit, as is usually the case. For this arrangement, 
solutions with an arbitrary contraction or expansion of the 
jet cannot exist, so that no effort need be made to avoid 
them. The solution for the general unsym metrical case, 
however, will show the jet velocity downstream at infinity 
to be different from the velocity upstream in the closed part. 
The possibility that, in the two-dimensional case, different 
pressures might be assumed on the two free surfaces still 
exists for this type of tunnel, but the resulting jet will have 
a constant curvature after leaving the neighborhood of the 
body. 

An upstream condition for the “infinitely long 1 ’ open 
tunnel and a correction to the results of reference 8. — In 
many discussions of the two-dimensional open tunnel, the 
set of images indicated in figure 5 (a) is used to satisfy the 
boundary condition that u~0, and the resulting flow shows 
an upflow in front and an equal downflow in back with no 
induced downwash at the wing itself. Actually', however, 
if the jet issues from a horizontal closed entrance — no matter 
how far upstream — it will remain essentially horizontal 
(because it is not subjected to any vertical force) until it 
reaches the wing. (See fig. 5 (b).) In order to eliminate the 
undesired upstream upwash, a uniform downwash should there- 
fore be added to the solution indicated in figure 5 (a). (Com- 
pare reference 7, p. 304.) Addition of this downwash does not 
affect the boundary conditions, since u is still zero at the 
boundary. This case is discussed quantitatively in part II, 


c 



(a) Deformation due to image system; no downwash. at the wing itself. 

(b) Undeformed upstream flow; downwash at wing is half of that at infinity. 
Fig cue 5. — Two-dimensional open tunnel of infinite length. 


where it is shown that the entrance-lip condition automati- 
cally provides the correct answer. 

Among the rectangular wind tunnels for which corrections 
were given in reference 8 is a type with closed sides but open 
top and bottom. The calculated corrections for approxi- 
mately square cross sections are approximately equal to those 
for the completely closed tunnel, a surprising result in view 
of the absence of any top or bottom constraint. The result 
is actually in error, as was discovered in an experimental 
effort to verify it (reference 9). In seeking to explain the 
errors the author of reference 9 pointed out that the image 
system used in reference 8 should have included an infinite 
row of vortices at infinity, and he showed how, by taking 
into account this row of vortices, the correct answer could 
be obtained. It could not be shown, however, that the 
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extent of this row of vortices is of a higher order , of infinity 
than is their distance from the origin, as is necessary if their 
effect is to be considered. The method of the preceding 
paragraph thus appears to he. much simpler and more rigor- 
ous in such cases than is a discussion of the image vortices 
at infinity. One simply observes that the image system of 
reference 8 provides an angle correction factor 8 of 0.25 for 
the flow far upstream of the wing, whereas 8 should be zero 
far upstream; a correction of —0.25 should therefore be 
added to all values of 8 computed by this image system for 
points within the tunnel. 

Summary of boundary conditions. — A basic physical 
characteristic of the flow is provided by the condition that 
the velocity be continuous at the entrance lip, which con- 
dition also helps to. provide uniqueness. The velocity on 
the free surface is not necessarily the velocity far upstream 
in the closed portion; in fact, for the two- d imensional case, 
it is even possible for the pressures on the two free surfaces 
to be different from each other. Equality of the velocities 
in the upstream and downstream closed portions has been 
recognized as an additional condition. Neglecting the upper 
portion of the closed exit may be desirable if the flow is so 
depressed that it does not make contact with the upper 
part of the exit. Neglecting the entire closed exit region 
may appreciably simplify the problem without introducing 
excessive inaccuracy if the region of interest is much closer 
to the entrance than to the exit. In general, adequate 
treatment of the exit (for large lift on the body in the tunnel) 
seems very unlikely. 

The discussion in the preceding sections has concerned 
mainly the physical flow conditions, and relatively little 
interpretation in terms of boundary conditions on the per- 
turbation potential has been given although such formal 
interpretation would appear a trivial task. The reason that 
this extension has not been made is that in a number of 
instances (as will appear subsequently) slight modifications 
of the basic viewpoint, leading to somewhat modified bound- 
ary conditions, are desirable for convenience of solution. 
Accordingly, the statements of the boundary conditions on 
the perturbation potentials will be given when the solutions 
are discussed. 

SUGGESTED ELECTRICAL ANALOGIES 
VELOCITY-POTENTIAL ANALOGIES 

Basic concepts of the analogies. — In the analogies to be 
discussed in the present section (none of which have yet been 
constructed), the perturbation velocity potential in the 
space within the wind tunnel is considered analogous to the 
electrical potential in a dilute electrolyte solution contained 
in a vessel of the same shape. An insulating material such 
as Bakelite, the conductivity of which is negligible compared 
with that of the solution, provides a boundary where the 

2 )$ 

normal potential gradient ^ is zero; and a metal, the con- 
ductivity of which is practically infinite relative to that of 


the solution, serves as a constant-potential boundary along 

which the longitudinal gradient ^ is zero. In such a setup, 

current is analogous to velocity except for a difference in 
sign (in the usual convention, current flows down a voltage 
gradient whereas air flow's up a velocity-potential gradient); 
in order to remove this difficulty, the sign convention for 
electrical potential is reversed in the following discussion. 

For greater clarity of exposition, the two-dimensional 
analogies are treated in detail, the three-dimensional 
analogies appearing as reasonably obvious extensions or 
modifications. It will be remembered, however, that any 
application will be found in three-dimensional problems 
inasmuch as most of the tv'o-dimensional problems can bo 
solved analytically. ♦ 

A two-dimensional vortex may be represented by two long 
metal plates separated by a thin insulator. (See fig. 6 (a) .) 
A flow corresponding to a vortex located at the edge of the 
plates is set up by applying a difference of potential across 



(b) 




. (a) Two-dimensional vortex showing current lines between the two plates. 

(b) Two-dimensional vortex in a perturbation field having a horizontal velocity component, 
(c) Three-dimensional element of lift. 

(d) Horseshoe vortex of finite span. 

Figure 6. — Velocity-potential analogies for two-dimensional and three-dimensional lifting 

elements. 
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the plates. If the perturbation, flow that results from, the 
presence of the vortex in the tunnel has a horizontal velocity 
component, this representation is no longer adequate because 
it requires the potential to be uniform along each plate. 
Rigor in this case would require that the plates be composed 
of a number of separate sections with each pair separately 
activated. (See fig. 6 (b) .) In this way it is possible to 
provide a potential difference between upper and lower 
surfaces that is everywhere equal to the desired circulation 
without requiring that the potential be uniform along the 
entire upper surface or lower surface. A horizontal velocity 
component normally occurs only when the lifting vortex is 
asymmetrically located in the tunnel. For simplification, 
only the simpler representation of figure 6 (a) is used in the 
remaining sketches. 

The element of lift in three-dimensional flow is the horse- 
shoe vortex of zero span, which is the same as a semi-infinite 
line of doublets. It may be represented by a pair of long 
narrow metal strips separated by an insulator. (See fig. 
6 (e).) As in the two-dimensional analogy, if the lifting 
element is asymmetrically located in the field, the strips 
must be made up of short pieces with each pair separately 
activated. The horseshoe vortex of finite span is repre- 
sented as in figure 6 (d), provided that there are no appre- 
ciable perturbation velocities in its plane. 

Evaluation of interference velocities. — The vertical veloc- 
ity component in the tunnel corresponds to the vertical 
voltage gradient in the electrolyte, which can be determined 
by measuring the voltage difference between a pair of short 
wire electrodes mounted one above the other a fixed distance 
apart. The tunnel interference at any point is found by 
measuring this voltage difference (relative to that across the 
two plates representing the vortex) first in the simulated 
tunnel and then in a large tank for which the boundary 
interference is either negligible or so small that it can be 
adequately computed by simple methods. Since the theo- 
retical flow field for the second case is known, the ratio of 
these two gradients, together with the distance from the pair 
of wires to the lifting vortex, should suffice to evaluate the 
boundary interference. The distance between the pair of 
wires need not be measured because only the ratio of the 
gradients is required. Similarly, the exact design and 
dimensions of the simulated lifting vortex are of no signifi- 
cance, provided that the gradients are determined at reason- 
able distances from it. In general, boundary interference 
at the vortex itself cannot be found directly by this method 
but may be determined by interpolation between or extra- 
polation from the neighboring points. 

Two-dimensional closed-open tunnel. — For simplification 
of the nomenclature, the open tunnel with closed upstream 
region but without a closed exit is designated the closed-open 
tunnel. The open tunnel with closed upstream and down- 
stream regions is designated the closed-open-closed tunnel. 

Figure 7 (a) illustrates the setup for a two-dimensional 
closed-open tunnel with a vortex on its center line. Shaded 
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(b) 


(b) Vortex on the center line. 

(b) Unsymmetricaliy located vortex. 

Ficttbe 7. — Velocity-potential analogies for the two-dimensional closed-open t unn el 


lines indicate msulatingboundaries ^ where -^=0^ andheavy 


unshaded lines indicate metal boundaries on which $ is 
constant. The upstream closed portion should be so long 
that the potential is essentially uniform at its upstream end; 
the length indicated on the figure should suffice. The open 
region should similarly be so long that the vertical flow 
between the vortex strips and the boundary strips no longer 
changes with distance downstream; again, the length indi- 
cated on the figure should suffice. From the condition of 

, . . . £>$ 
velocity continuity at the entrance lips and the fact that ^ 


is zero along the free boundaries, it follows that ^ must be 

zero at tbe edges of the two closed boundaries. Tbe poten- 
tials on the two free boundaries must, therefore, be adjusted 
until the difference between the potential of each and the 
potential of a thin feeler electrode just upstream of its edge 
is zero. For the symmetrical condition shown, the single 

variable voltage source indicated will provide zero ^ 


at both edges simultaneously. 

Figure 7 (b) illustrates the setup for the two-dimensional 
closed-open tunnel with the vortex in an off-center position. 
A single variable voltage source across the two free bound- 
aries is now no longer capable of simultaneously satisfying 
the continuity condition at both edges, so that an additional 
variable voltage source and an upstream electrode are 
required. Tbe current in the closed part of the t unn el flowing 
into this upstream electrode corresponds to an upstream 
perturbation velocity. This upstream perturbation velocity 
constitutes tbe previously mentioned difference between the 
velocity far upstream in the closed part and the velocity on 
the free surface. The concept here is slightly at variance with 
previous discussion, which considered a perturbation velocity 
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along the free surface, with the far upstream velocity appear- 
ing as the undisturbed velocity V. As the analogy is set up, 
however, no perturbation velocity may appear along the free 
surfaces because they are at constant potential; hence, the 
total velocity on the free surfaces must be considered as the 
undisturbed velocity U and any difference between this 
velocity and the velocity far upstream appears as an up- 
stream perturbation velocity. As appears in part II, this 
viewpoint is also found convenient in the analytical solution 
of these problems. 

Two-dimensional closed-open-closed tunnel. — The setup 
for the two-dimensional closed-open-closed tunnel with a 
vortex on the center line (fig. 8 (a)) is an obvious modification 
of the corresponding setup for the closed-open tunnel. The 
same would be true for the off-center vortex except for the 
necessity of satisfying the condition that the velocities in 
the upstream and downstream closed regions be'equal. Thus, 
the setup of figure 8 (b) provides an upstream perturbation 
velocity but no downstream perturbation velocity and can- 
not, therefore, solve the problem completely. An additional 
flow, found by the setup of figure 8 (c), must be included. 


/ ///////////////////////, 





(a) Vortex on the center line. 

(b) Unsymmetrically located vortex (incomplete representation) . 

(o) Expanding or contracting Jet. 

Figure 8. — Velocity-potential analogies for the two-dimensional dosed-open-closed tunnel. 


An electrode is here located at both the upstream and down- 
stream ends, and the potentials relative to the free boundary 
are so adjusted that the entrance-lip condition is satisfied. 
It is apparent that in order to satisfy this condition the 
downstream current flow will be much greater than the 
upstream current flow; that is, the downstream perturbation 
velocity for a contracting or expanding jet is much greater 
than the upstream perturbation velocity. Because of this 
difference, a suitable amount of the flow of figure 8 (c) may 
be added to that of figure 8 (b) to produce equal upstream 
and downstream perturbation velocities. 

Displacement of the free surfaces. — The current density 
normal to the surface of a metal plate representing a free 


surface is proportional to ^ and corresponds to the local 
vertical perturbation velocity. The total vertical displace- 
ment at a point on the free surface is then given byj*^ dx 

integrated from the entrance lip to the point. In particular, 
the integral along the entire lower free surface_of a closed- 
open-closed tunnel represents the displacement at the exit 
lip and it may be measured by means of an ammeter in the 
line that goes to the lower metal plate. 

If the pressure on the lower free surface can adjust itself 
so that the displacement at the downstream end is zero, the 
longitudinal perturbation velocity at the lower surface will 
be different from that at the upper surface. If the perturba- 
tion velocity on the upper surface is taken as zero, that on 
the lower surface will be negative, so that the potential on 
the lower surface must drop uniformly from entrance to exit. 
Such a variation could be accomplished if the lower surface 
were represented by a number of short metal strips instead 
of a single plater - 

Three-dimensional closed-open and closed-open-closed 
tunnels. — The analogies for the three-dimensional tunnels 
are obvious modifications of those for the two-dimensional 
tunnels. The free boundary may not be simulated by a single 
metal cylinder because, as was previously noted, different 
elements of the free boundary do not have the same potential, 
although the potential is constant along each element. The 
free boundary must thus be simulated by a number of longi- 
tudinal metal strips, insulated from each other, with a feeler 
electrode immediately in front of each. When the lifting 
element lies in the horizontal plane of symmetry, the en- 
trance-lip condition may be satisfied without an additional 
electrode. If the lifting element is not in the horizontal 
plane of symmetry, an upstream electrode will be needed, 
with the potential of each strip adjusted relative to this 
electrode. For the unsymmetrical closed-open-closed anal- 
ogy, the requirements that upstream and downstream 
velocities be equal necessitates further measurements with a 
setup corresponding to that of figure 8 (c). 
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ACCELERATION-POTENTIAL ANALOGIES 

Basic concepts of tlie analogies. — The pressure has the 
properties of a potential — designated acceleration potential — 
in a field consisting of a small perturbation flow superposed 
on a uniform stream. If the pressure in the undisturbed 
stream is taken as zero, then the perturbation velocities are 
related to the pressure by the following equations: 
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where 

p density 

t time 

p pressure 

Since, by the first equation, u is proportional to p, it is simpler 
merely to consider the perturbation velocity u itself as the 
potential, with o and w given by the following equations: 


and 





The necessity of performing an integration in order to de- 
termine v or w is a basic disadvantage of the acceleration- 
potential analogy compared with the velocity-potential 
analogy in which o and w are measurable directly. 

In the analogies to be discussed in the present section, the 
perturbation velocity u is considered analogous to the elec- 
trical potential in a dilute electrolyte solution. A metal 
serves as a boundary along which u is constant, and the local 

intensity of current flowing into it gives an insulator 

serves as a boundary where ^ is zero. The lifting element 

in either two or three dimensions is represented by a pair of 
short metal plates separated by an insulator; when the upper 
plate is maintained at a higher potential than the lower plate, 
the arrangement represents a thin airfoil with suction (large 
u) on its upper surface and pressure (small or negative u) on 
its lower surface. The current at each of the two plates 
should be the same in order that the slope of the airfoil sur- 
face (proportional to v) be the same on both upper and lower 
sides. In order always to satisfy this condition the voltage 
source activating the lifting element should not be tapped 
to any other electrode in the field. 


is 


Two-dimensional closed-open tunnel. — The setup for the 
two-dimensional closed-open tunnel with the lifting element 
on the center line is shown in figure 9(a). The walls of the 
upstream closed region are represented by insulators, which 

establish that ^=0 at every point; hence, the condition 

that »= J ^dx = 0 at every point on tbeclosed boundary! 

satisfied. The two free boundaries are represented by metal, 
and electrically connecting them, as shown, satisfies the 
further condition that they have the same potential (the same 
u ) . The flow of current into the lower boundary then equals 
the flow of current out of the upper boundary, so that the 
ultimate downstream value of v will be the same on both up- 
per and lower boundaries, as is desired. In fact, for the 
symmetrical case illustrated, the value of e will be the same 
at all pairs of opposite points on the two free boundaries; 
so the boundaries will be everywhere parallel. No special 
attention need be paid to the entrance lips; the entrance-lip 
condition is automatically satisfied since the potential u is 
continuous at these points (although the potential gradients 
at these points are infinite) . 

For the off-center position of the lifting element (fig. 9 (b)) 
no modification of the circuits is needed. The difference 
between the potential in the upstream closed region and the 
potential of the free boundaries, which is the upstream per- 
turbation velocity u, is measured with the aid of the probe P. 
As in the symmetrical case, the ultimate downstream value 
of v wrill be the same for both the upper and the lower bound- 
aries; however, it is no longer true that the two boundaries 
will be everywhere parallel, and the ultimate width of the 
jet will be different from the width of the closed part. 

Two-dimensional closed-open-closed tunnel. — The value 
of v at the downstream end of the free boundary is given by 
the total flow of current into the metal plate and, in general, 
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(a.) 



(b) 


(a) Lifting element on the center line. 

(b) Unsvmmetrical location of the lifting element. 

Figure 9. — Acceleration-potential analogies for the two-dimensional closed-open tunnel. 
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is not zero. At the lip of the closed exit, however, v must 
be suddenly reduced to zero in order for the flow to follow 
the solid boundary; hence, a short electrode must be added 
at the exit lip, and as much current must be forced out of it 
as flows into the long electrode that represents the open 

boundary — that is, the integral of ^ dx along the free bound- 
ary must be canceled at the exit lip. The setup (fig. 10 (a)) 
therefore shows a voltage source to supply this current and 
means for measuring and equalizing the current flow into 
adjacent electrodes. If these additional short electrodes 
are omitted, the setup will correspond to a tunnel the exit 
section of which has been alined with the deflected jet (fig. 

10 (b)) because the condition that ^=0 on the closed exit 

boundary would merely permit v to remain at the value it 
had at the end of the free boundary. 
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(a) 




(b) 



(c) 



(d) 


(a) Lifting element on the center line. 

(b) Tunnel arrangement that corresponds to omitting the additional short strips. 

(c) Expanding or contracting jet. 

Cd) Curving jet. 

Figtjbe 10.— Acceleration-potential analogies for the two-dimensional closed-open-closed 

tunnel. 


For the off-center position of the lifting surface, a similar 
setup is used and, as before, probes in the regions far up- 
stream and far downstream are used to determine the poten- 
tial u in these regions relative to the potential of the free 
boundary. Since these potentials far upstream and far 
downstream will not be equal, an additional perturbation 
field must be provided such that the sum of the two fields 
will have the same potential in the two regions. This 
additional perturbation field, which corresponds to a con- 
traction or expansion of the jet, is provided by the setup 
shown in figure 10 (c) . It is clear from this figure that the 
downstream perturbation potential is much greater than the 
upstream perturbation potential; this result corresponds to 
that indicated in the velocity-potential analogy. 

The condition in which the pressure on the lower free 
surface is higher than that on the upper free surface is easily 
represented by applying a voltage difference between the 
two surfaces. (See fig. 10 (d).) The corresponding dis- 
placement of the lower surface, however, is not so readily 

obtained. The vertical velocity at every point is dx, 
so that the displacement at each point is //!?**■ In 

order to accomplish this integration ^ must be determined 

at points along the boundary, perhaps by breaking the long 
plate into a number of short pieces and determining the 
current flowing into each. 

Three-dimensional closed-open and closed-open-closed 
tunnels. — The analogies for the three-dimensional tunnels 
are again obvious modifications of those for the two-dimen- 
sional tunnels. For the closed-open analogy, the free bound- 
ary may be represented by a single cylinder of metal (fig. 
11 (a)). For the closed-open-closecl analogy, the free 
boundary must be represented by a number of separate 
strips (fig. 11 (b)) in order that the total current into each 
strip may be measured and an equal current forced out of the 
short strip immediately behind it. Contraction or expansion 
of the jet is represented as in the two-dimensional case, but 
the setup that would correspond to different pressures along 
different strips seems to have no practical significance in the 
three-dimensional case. 

Some form of this acceleration-potential analogy is prob- 
ably the most convenient for solving problems similar to 
that of the helicopter in the Langley full-scale tunnel. 
Simply neglecting the exit, as with a closed-open tunnel, 
permits the free surface to be represented by a single sheet 
of metal and eliminates any measurements of current flow to 
or from the surface. Improved accuracy should be attain- 
able by cutting the sheet into two parts with two short strips 
at the rear. (See fig. 11 (e).) The need for many strips 
seems unlikely, at least in view of the previously mentioned 
uncertain definition of the physical flow in the region of the 
exit. 

Correspondence between velocity-potential and accelera- 
tion-potential analogies. — As has already been indicated, 
the acceleration potential is identical with the ^-component 
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the closed-open-closed tunnel. 

of the perturbation Telocity and is hence merely the 
E-derivative of the perturbation-velocity potential. It is of 
interest to point out the related fact that the acceleration- 
potential analogies are, in a sense, the z-derivatives of the 
velocity-potential analogies. For example (see fig. 12), 

(1) For the velocity-potential analogy, an infinitely long 
double layer represented a lifting element located at its 
forward edge. The. difference between two such double 
layers, of which one is shifted slightly relative to the other, 
is merely the short double layer that was used in the 
acceleration-potential analogy. 

(2) For the velocity-potential analogy, the free boundary 
co nsis ted of constant-potential strips on which the potentials 
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(b) 

(a) Representation of lifting element. 

(b) Representation, of open boundary in closed-open-closed, tunnel. 

Figure 12. — Acceleration-potential analogies as the difference between two velocity-potential 
anal ogies slightly shifted relative to each other. 


were so adjusted that the gradient was zero at the leading 
edge. If each strip is now shifted and subtracted, there 
remains a long strip, with a short strip at the front and back. 
Since, in the velocity-potential analogy, the gradient was 
zero at the entrance lip, the short strip at the front may be 
neglected. The remainder corresponds to the arrangement 
used in the acceleration-potential analogy, and the fact that 
the total current after the subtraction must be zero corre- 
sponds to the fact that the total current out of the short 
strip must be made equal to the total current into the long: 
strip. 

(3) When the lifting element was off-center, the velocity- 
potential analogies required electrodes upstream and down- 
stream with uniform current flow along the upstream and 
downstream closed regions. That the subtraction eliminates 
these current flows corresponds to the fact that no upstream 
or downstream electrodes are used in the acceleration- 
potential analogies. 

TECHNICAL DIFFICULTIES 

It should be pointed out that the analogies herein described 
may be rather unwieldy, experimentally. Even for the 
simplest types of analogies, the literature indicates con- 
siderable uncertainty as to the most satisfactory electrolyte 
and electrode materials, and appreciable difficulty in balanc- 
ing capacitances (alternating current is generally used in 
analogies in order to minimize polarization at the electrodes). 
In the present analogies, the need for separate current sources 
that are exactly in phase and the large capacitances that 
will certainly characterize the vortex and the open-boundary 
representations should greatly complicate the technique. 
Perhaps the use of direct current instead of alternating 
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current, with nonpolarizing electrodes (as platinized plati- 
num), would be a more practical approach in this respect. 
Simultaneously satisfying the entrance-lip condition at a 
number of points around the inlet (or satisfying the cor- 
responding exit condition for the acceleration-potential 
analogies) may also turn out to be very difficult. 

RfiSUMfi OF PART I 

The most significant points of the preceding discussion of 
open wind tunnels and their electrical analogies are as 
follows: 

1. Continuity of velocity at the lip of the entrance cone is 
a basic characteristic of the flow in an open wind tunnel. 

2. Equality of the velocities in the upstream and down- 
stream closed regions is a further condition on the tunnel 
flow if extraneous longitudinal pressure gradients are to be 
avoided. 

3. The velocity on the free surface need not be the same 
as the velocity in the closed upstream region. In general, 
the two velocities are the same only when the lifting element 
lies in the plane of symmetry of the tunnel. 

4. For the two-dimensional open tunnel the velocities on 
the two free surfaces need not be equal. If the space below 
the lower free surface is closed off, the pressure on the lower 
free surface will adjust itself so that the displacement at 
the exit lip is zero. 

5. Considerable uncertainty exists with regard to condi- 
tions at the exit or the mathematical equivalents of these 
conditions. Correspondingly, certain compromises in com- 
plying with the idealized downstream conditions may be 
justified in a determination of boundary interference. 

6. In any analysis that neglects the closed entrance and 
exit regions, the condition of zero upstream induced flow 
must be retained. 

7. Electrical analogies of either the velocity-potential or 
the acceleration-potential type may be devised to correspond 
to most of the problems discussed. 

8. In electrical analogies that represent velocity potential 
by electrical potential, the condition of continuity at the 
entrance lip appears troublesome, especially for three- 
dimensional tunnels- however, the exit conditions are easily 
represented. 

9. In electrical analogies that represent acceleration 
potential by electrical potential, the entrance-lip condition 
is automatically satisfied but fulfillment of exit conditions is 
troublesome. Rough approximation of the exit conditions 
may, however, be adequate for many purposes. 

10. Acceleration-potential analogies are experimentally 
simpler than velocity-potential analogies. 


II. TWO-DIMENSIONAL TUNNELS 

In part II, boundary-induced velocities in two-dimensional 
open tunnels are'derived ■with special reference to the effects 
of the closed entrance and exit regions. Of special interest 
is the fairly simple mathematical procedure introduced for 
solving flow problems of this type. The cases treated are: 

(1) Tunnel with a closed entrance (upstream) region but 
without a closed exit region 

(2) Tunnel with a closed entrance region but with only one 
exit lip (corresponding to a condition in which the downward 
deflection of the flow is so large that the flow makes contact 
only with the lower exit lip) 

(3) Tunnel with closed entrance and exit regions 

(4) Same as case 3, but with different pressures on the 
two free surfaces 

Numerical results are given for all cases. 

SYMBOLS AND DIMENSIONS 

Each tunnel is idealized as a strip of uniform height k, 
having a stream velocity V, and containing a point vortex 
of strength r. For simplification of the present develop- 
ment, lengths and velocities will be made nondimension al 
by dividing by h and V, respectively, and the vortex strength 
will be made nondimensional by dividing by hV. Essen- 
tially, then, the solutions will be developed for a vortex of 

p 

strength ?'~Ty * n a tunnel of unit height; and in the follow- 
ing list of symbols the lengths in the complex planes arc 
in terms of h, and the complex velocities are in terms of V: 

h tunnel height 

V tunnel velocity 

f complex variable of physical plane (J-f-iij) 

location of vortex in ("-plane 

s complex variable of transformed plane (a: -HjO 

Zi location of vortex in s-planc 

2 complex velocity in physical plane (u—iv) 

Q complex velocity in transformed plane (u — iv) 

A, B, C, M,N real constants 
r 

r 

2<(ri) 
a 

Wi • 


vortex strength 

nondimensional vortex strength ) 

induced complex velocity at f when vortex 
is at fx 

induced complex velocity at f, 
abscissa of exit lip in transformed space 
a complex velocity in the form of an elliptic 
integral of the first kind 
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a comples velocity in the form of an elliptic 
integral of the second kind 
a comples velocity 

l variable of integration 

G function defined by equation (2) 

K,K' complete elliptic integrals of the first kind, 

with modulus 1/a 

E, E' when not followed by parenthesis, complete 

elliptic integrals of the second kind, with 
modulus 1/a; with upper limit indicated 
in parentheses, incomplete elliptic integrals 
of the second kind, with modulus 1/a 

F, F r incomplete elliptic integrals of the first kind, 

with modulus 1 fa, and with upper limit 
indicated in parentheses 
R-P- real part 

1- P • imaginary part 

c airfoil chord 

Ci airfoil lift coefficient 

e tunnel-induced downwash angle, radians 

horizontal perturbation velocity at free 
boundary 

BOUNDARY CONDITIONS 

The two-dimensional tunnels discussed are considered to 
have their fixed and free boundaries parallel to the real axis 
with the main tunnel flow from left to right. The physical 
plane (in which lengths and velocities have been made non- 
dimensional as Just described) will be designated the f-plane, 
with the complex perturbation velocity u~iv, or q(f), subject 
to the following conditions: 

(1) On each fixed (or closed) boundary, I.P. q(f)s—y=0. 

(2) On each free (or open) boundary, R.P. q/^)=u=0 or a 
constant. 

(3) At each lip of the closed entrance section, g(f) is 
continuous. 

CASE I— CLOSED-OPEN TUNNEL 

Total perturbation velocity. — By the transformation 

z=e*s ( 1 ) 

the tunnel in the f-plane, represented by an infinitely long 
strip of unit height, is transformed to the upper half of the 

2 - plane. The correspondence between points is shown in 
figure 13. 


•*i 

Figure 13. — Physical and transformed spaces for bvo-dimensioaal closed-open 
tunnel of unit height. z=«*r. 

The complex velocity (rather than the more usual complex 
potential) is considered to be retained in the transformation, 
and the problem is thus to find a fimction Q{z), where 

u— i» = Q( s)=g(f) 

such that 

(1) On the closed sections of the boundary, that is, for z 
real and [z[<l, 

I.P.<2(z)=0 

(2) On the open sections of the boundary, that is, for z real 
and |z|>l, 

R.P.<2(2)=0 

(3) Forz=±l, 

Q(z )= 0 

(4) Q(s) is finite at infinity. 

Consider the complex velocity 6{z) corresponding to a 
vortex at z x and its reflection at z t : 

(2 > 

This function, which is of order 1/z 2 at infinity, satisfies 
conditions (1) and (4) but not (2) and (3). Functions of the 
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form Yl — z 2 , gVl — g 2 , . . . satisfy conditions ( 1), (2) , and 
(3). At infinity, Yl — g 2 is of order g and z-\ll —z 2 is of 
order g 2 ; therefore, either of the products G^l—z^Gz Yl— g 2 , 
or a linear combination of the two satisfies the four condi- 
tions and has a pole of the first order at z v These facts 
suggest that the desired velocity function is of the form 

e<°)=<^-^S>+*>Vl=? <3) 

where A and B are as yet undetermined real constants. 

The values of A and B are to be determined such that the 
pole at g t represents a vortex of strength V in the f-plane. 
Thus 


r'= j> 2 (f) dt = <j> Q{e) J dz =i ^ Q(z) ~ 

where the integral is taken about the point gi. By Cauchy’s 
integral formula 

T'=-2{A+Bz x ) (4) 

#1 


where the last equality follows from equation (4). The 
remainder of equation (5), after 7r(f — f j) is replaced by 

log is 


i 

(A +Bz) Yl -2 s log •~-+W- (z — 2 i) 
i g-L „f. - 

(g— gOlogj- 

This expression is of the form ^ for g = g,. Differentiating 

numerator and denominator, according to L’Hdspital’s rule, 
still leaves both equal to zero at z=z t (that the derivative of 
the denominator is zero at z—zi is obvious; that the deriva- 
tive of the numerator is also zero at g^z, can be verified with 
the aid of equation (4)). A second differentiation yields the 
following expression for the limit as z approaches z t : 


{A-\-Bz)-\!T—z i , r 7 


Z — Zl 


2 log— 

2l_ 


The values of r 7 and z : are known, so that this complex 
equation can be solved for the two real constants A and B. 
Substituting these values in equation (3) will thus give the 
desired complex velocity function. 

Tunnel-interference velocity. — The tunnel-interference ve- 
locity is defined as the difference between the total perturba- 
tion velocity 2 (f) due to the presence of the vortex in the tun- 
nel and the velocity due to a vortex in an unbounded medium; 
that is, the tunnel-interference velocity 2f(f,fi) is 


2 ,<r, ro=< G^-jYr,) *+*+*=*+*£& 

(5) 

If the vortex is on the axis of the tunnel (that is, if 
then from equation (1), Zi=iy u and 
equation (4) gives 

A=0 
r 7 

B= ± 

2Yl+lh 2 



If the point of evaluation is also on the axis (that is, 
£ = £+ whence z=iy\ then 


!< 0+5’ &+|) — ir ' pcr^Yl 


Iy 2 + 1 

, ir 7 

\lyi 2 + 1 

' 2 


(7) 


Thus, if the vortex is on the axis, the interference velocity at 
all points on the axis has only a vertical component. 

The interference velocity at the vortex itself is the limit of 
expression (5) as z approaches gi. The term containing g— 1 X 
offers no difficulties and its limit is readily evaluated: 


lim— — — ( A+Bz ) YI^ 2 = ~~ (A+BzJ 

g — 2 1 4 c y 1 


— (A+5gi) Yh— g? 

i 


2(A+i?z,)g t 2 

Vi— g t 2 

2Zi 


2Bzi-\h — z l i 


This fraction can be greatly simplified by use of equation (4), 
and the result, added to the previously derived limit, gives 
the desired correction at the vortex: 


Iha 2<(fj fi)=2i(fi)=i - 

r-fi L 


r . 


r 7 


4 1 2(l-g 1 2 ) 


-B^l—zi 2 - 


r 7 g t 1 


4'iyJ 

(S> 


For the special case in which the vortex is on the tunnel 
axis ( z=iy ), this expression reduces to the following form 
(after substituting for B from equation (6)): 



v 1 2 

1+2/1 2 


(9> 


Upstream perturbation velocity. — If the vortex is not on 
the tunnel axis, A will not be zero (compare equation (6)). 
The tunnel interference velocity far upstream in the closed 
part of the tunnel is found by putting £= — <» and g=0 in 
equation (5), which then reduces to 


2 a-*, (A -£)■ 


2Ay { 


Zl 


which is real. For this unsymmetrical case, therefore, a finite 
longitudinal perturbation velocity is found far upstream in 
the closed part. As was pointed out in part I, such results 
appear because the problem was set up so that the longitudinal 
perturbation velocity on the open boundary is zero. If 
the velocity far upstream in the closed part is to be taken 
as the base, the result means merely that the longitudinal 
velocity on the open boundary exceeds this base velocity by 

• j —p 1 and that there is a corresponding difference in pressure 
between the closed part and the space surrounding the jet. 
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limiting case of completely open tunnel. — With increasing 
distance of the vortex from the closed entrance (that is, 

ir' 

with increasing iji), expression ( 9 ) approaches — 5— As was 

pointed out in part I, the image system normalLy used to 
satisfy the boundary condition on an infinitely long, open, 
two-dimensional tunnel produces no induced flow at the 
vortex itself, and only after introduction of the additional 
condition that the upstream flow be horizontal is this value 
■iT' 

of — for the induced-flow correction obtained. In the 

present development, however, it is seen that the condition of 
continuity at the entrance lips automatically takes care 
of this condition on the upstream flow direction, even when 
the entrance and the vortex are infinitely far apart. 

No further discussion of the completely open t unn el will be 
given here inasmuch as this case has been adequately 
treated by the method of images. (See reference 6, p. 302 .) 

CASE 2— TUNNEL WITH ONE FIXED EXIT BOUNDARY 

Perturbation, velocity. — As before, the transformation 
2=e ir transforms the tunnel, considered as an infinite strip 
of unit height, into the upper half of the s-plane. The 
correspondence between points is shown in figure 14 . The 
conditions on the complex velocity 0 (z) are: 

(1) On the real axis, I.P.O(2)=0 for — 1<2<1 and 
for 2>a. 

(2) On the real axis, R.P.O(s)=0 for 2<— 1 and for 
!<-<&• 

( 3 ) For2=±l, (?(z)= 0 . 

( 4 ) Q{z) is finite everywhere in the upper half-plane except 
at z=a and at 2=2j. As noted in part I, 0 {z) will be infinite 
at z=a. 
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Figure 14.— Physical and transformed spaces for two-dimensional tunnel of unit height 
with one exit boundary. z=& r i. 


The function G(z), given by equation (2), will again be 
used as a factor that is real along the entire real axis and 
has the desired type of singularity at zj. The functions 


jl—2 2 ll — z 2 

-t/ _ - and z-\j change from pure real to pure imagin- 

Y cl z \ a z 

ary as z passes through ±1 or through a, and, furthermore, 

are zero at s=± 1 and infinite at z=a. At infinity they 

are of order 2 1/2 and z 3/2 , respectively. Therefore, as before, 

L — 2 2 , nf s ll—z 2 

—and G[z)z-t 

2 V a— 2 

satisfies the preceding conditions and has a pole of the first 

order at Zi. Q(z) is therefore of the form: 


a 


jl—z 2 ll—z 

linear combination of G(z )-\ / — — -and G{z)z-\ 

' \ a—z V 0.—1 


< 10 > 

The real constants A and B are determined by the same 
condition as before, which gives 


r' 


2{A-{-Bz\) 1 1 — z 2 
2i V a—z t 


( 11 ) 


Tunnel-interference velocity. — The interference velocity 
2i (r,ri) is the difference between the perturbation velocity 
and the velocity due to a vortex located at the same point in 
an unbounded medium: 


(i=iri = %) ' A + Bs > VS 


iT' 


z ^2x(r-r.) 


( 12 ) 


This expression does not simplify appreciably if the vortex 
is located on the axis, and the interference velocity at a point 
on the axis due to a vortex on the axis is not normal to the 
axis. 

The interference velocity at the vortex itself is the limit 
of expression (12) as 2 approaches z x . Proceeding as in the 
preceding case gives 



T'a n / I , r 'z t 
Ua-zr \ a — 2i _r 4i]/ 1 _ 


( 13 ) 


It may be shown with the aid of equation (11) that this equa- 
tion reduces to that for the closed-open tunnel as a goes to 
infinity. 


CASE 3— CLOSED-OPEN-CLOSED TUNNEL 

Perturbation velocity. — The transformation z=« xf trans- 
forms the tunnel space into the upper half of the 2-plane with 
correspondence between points as indicated in figure 15 . 
The boundary conditions on the complex velocity Q{z) are: 

(1) On the real axis, I.P.(?(2)=0 for |z|<l and for [z[]>a. 

(2) On the real axis, R.P.< 2 (s )=0 for l<CM<C a - 

( 3 ) For z—±l,Q(z)= 0 . 

( 4 ) Q(z) is finite everywhere in the upper half-plane except 
at 2= ±a and at z=z x . 

(5) <?(0) = <?(°°), (this equation corresponding to the 
condition noted in part I that the perturbation velocities 
in the upstream and downstream closed regions he the same). 
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B' 


B' 


$ -plane 
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r 

B 




z-p/one 


i C' D' 




iy 

B' a\ 
-1 


' A B rr D_ 

~Z n 


X +1 


Figure IS.— Physical and transformed spaces for symmetrical two-dimensional elosed-open- 
ciosed tunnel of unit height. z=« T I. 

The function 6(z) given by equation (2) is again used as 

VI •• 2 ^ j \ g ® 

a 2 ~ z 2 ‘ Z ~\ a 2_ g P an< ^ 

V \ g2 

— 1 satisfy conditions (1), (2), and (3) and are of 

a z 

orders 1, s, and z 2 at infinity, respectively. By the same 
reasoning as before, 


< 14 > 

The condition that the pole of the first order at Z\ represents 
a vortex of strength r' is 


r=- 


2(A+Bzi+Cz 


z 1 


i2 I lull 
M a 2 - 01 * 


(15) 


Condition (5) is satisfied by equating the two forms of 
equation (14) for z equal to zero and equal to infinity. Thus 


Q(0)=i 

V *t a 

... — 2yi -A 
Nil 2 a 


. lim Q(z)=\unj——^-=r(A+Bz+Cz 1 ) 

Z— >C3 Z— >a> ^1/ \ & & 


=—2y l G 

whence, by condition (5) 

— 2yi A_ 


Nil 2 a 


-%ViC 


or 


<7= 


a|si| 2 


(16) 


The complex equation (15) and the real equation (16) suffice 
for evaluating the three real constants A, B, and C in 
equation (14). 

Tunnel interference velocity. — The tunnel interference 
velocity is 


»r 


2*tt-fi) 


(17) 


If the vortex is on the tunnel axis, that is, if Zi=iy u equa- 
tion (15) gives 


B ~~i l+yl 


< 7=4 

Vi 


(18) 


Comparing this last equation with equation (16), which 
reduces to the following form for Zi=iyi 

„ «Vi 

shows tliat for this symmetrical case 

A=C^ 0 

since, as is clear from figure 1 5, a 1 . The interference veloc- 
ity at a point on the axis (z—iy) due to a vortex on the axis 
is thus 

- A. , i i A ir'iM l 1 +y 2 « 2 +i/i 2 , ir' , . 

which is normal to the axis. 

The interference velocity at the vortex itself is the limit 
of expression (17) as z approaches 2 i. Proceeding as before 
gives 

.["T'.rw a 2 — 1 

2 ( a *-g l «)( 1 -g l *) + 

< s + 2 WiS? + £§ rO (20) 

. For the case in which the vortex is on the axis (z\=iyl 
the normal velocity at the free boundary (z=x, where 
l<a:<a) is given by 


Q{x)=i(—. -44 BxJlf- 

\x—iyi x+iyj V a 2 — 


— 

X 2 


or 


Q(x)-- 


2iy± 


x 2 +y i 2 




( 21 ) 
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CASE 4— CLOSED-OPEN-CLOSED TUNNEL WITH UNEQUAL 
PRESSURES ON THE FREE SURFACES 

Boundary conditions. — As indicated in part I, the two- 
dimensional closed-open-closed tunnel may develop unequal 
pressures on the two free surfaces if a closed space exists 
below the lower free surface. Within the limits of the 
present linear theory, this pressure difference corresponds to 
superposing on the flow discussed in the preceding section an 
additional perturbation velocity field Q(z) that 

(1) Has no singularities within the t unn el 

(2) Satisfies the condition of continuity at the inlet lips 

(3) Has a horizontal component equal to, say, -f 1 on the 
upper free boundary and —1 on the lower free boundary 

(4) Has no vertical component on the closed boundaries 

(5) Is zero at infini ty upstream and downstream 
Rewritten as conditions on the complex velocity Q(z) in 
the z-plane (fig. 15), these conditions become: 

(1) Q(z) has no singularities in the upper half of the 
z-plane. 

(2) <2(z)=0 at z=±l. 

(3) On the real axis, Q(z)=—1 for l<z<a, and 
Q(z) = + 1 for — a<z<— 1. 

(4) On the real axis, I.P. Q{z)— 0 for |z|<l and for 
|z[>a. 

(5) Q(z)=0 for z=0 and for 2 = 00 . 

Outline of method. — Consider the following two functions 
of z: 

m = r z dz 

Wl Jo 


i\f x-plone 

i 

QL D' r r B' A'U B r r D C 

~ -« Jr 1c 7i * * a =* 


D’ C'C D 


B' A'\A B 


Wife) 



velocity function for the closed-open-closed tunnel with 
unequal pressures is thus of the form: 

Q(z) — Mw^iz) +Nw 2 {z) +w 3 {z) 



They can be considered as complex velocities having the fol- 
lowing properties along the real axis (compare reference 10) : 
Wi is real between 1 and — 1 ; between 1 and a or between 
— 1 and —a, its real part is constant but an imaginary part 
is introduced; beyond a or —a, the imaginary part 
is constant while the real part approaches zero; 
R.P.u- 1 (2) = -R.P.w l (-2); I.P.«j 1 (z)=I.P.to I (— 2 ). 

v: 2 has the same properties as w t except that beyond a 
and —a its real part approaches <= and — <= , respectively. 

Maps of the two functions are shown in figure 16. It 
should obviously be possible to find a linear combination of 
these two functions, Mvh+Nwz, such that for l< 2 <a 

R . P . (Mwi+Nvjz) = — 1 

for — a< 2 <— 1 

R. P. (Mwi-\-Nw 2 ) = + 1 
and beyond a or —a 

I. P. (Afw 1 -i-iV r W2)=0 

A simple additional function w 3 to be discussed subsequently 
is needed to satisfy the condition at infinity. The desired 


The constants M and A r are derived in the two following 
sections. 

Evaluation of integrals. — In the following development, 
the modulus of all the elliptic integrals is 1/a; the modulus 
t will therefore not be indicated in the symbols. In the 
designations for the incomplete elliptic integrals, E, E', F, 
and F', the terms in parentheses are the upper limits of 
integration. Then 

-d)—(«)=r va K 

Wo(l) =R.P .w 2 {a ) = J y ° dz=aE 

* 

I.P- Wl(a)= Ji V(l-e}(a f -z5 

which by the substitution z 2 =a 2 — (a 2 — 1)Z 2 reduces to 



I.P.U' 2 (a)=J^/— ^dz 
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which by the substitution z 2 =a?—(a 2 —l)l 2 reduces to 


-r 


a 2 - 1 


h n 

a 2 ■ n <ti 


—ia 



1 - 2 ^= 1 1 > 

a 2 

l-Z 2 


■ dl =iaK'—iaE' 


Solution of simultaneous equations for .M and iV. —With the aid of the four formulas just derived, the two previously 
mentioned equations in M and N may be written 

— K+NaE=-l 
a 

— K'+NaK'-NaE'=0 
a 

which are easily solved simultaneously for M and N. By introducing the following relation between the complete elliptic 
integrals (reference 11, p. 520) 

EK'-KK'+KE'^ 

the expressions for M and N are finally obtained in the following forms: 

M=— (K'-E>) 

7 r 


N= 


2 K' 


<Ltc 


Value of Mw,, + Nw 2 at infinity. — The constants M and N have been determined so that I.P . (Mwi +Nw 2 ) — 0 at infinity; 
furthermore, R.P.Mwi— 0 at infinity, as is clear from figure 16. Therefore, the value of Mwi+Nw 2 at infinity is merely 
R.P.A %2 at infinity. It is necessary to investigate this limit before choosing the form of w 3 , because, as was previously 
noted, the purpose of w 3 is to provide Q(z )= 0 at infinity. The limit may be written 

The first term is simply NaE. In order to evaluate the second term, substitute z=y - 

—Na 2 


■/: 


- 1 2 1 2 
i -i 2 


dl 


V(1 l 2 


— —ATn2 f 1 — -Jf-WnZ f 1 — 

Jo vcr=w=ly Jo i 2 Va-iW- 


1 2 VU -l*)(a 2 -l 2 ) 


The first term is —NaK. In order to evaluate the second term, it is noted (reference 12) that 
d V(l-^(« 2 -j 2 )_ _ Z*-o> _ -a 2 V 12 -<^ 2 
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Transposing terms in this equation gives 

a 2 a 2 ja^-P 

ZVCl-PKa 1 -^ V'(l— ^(a 2 — Z 2 ) V l-? 2 

rf V(l-^(a 2 -/ 2 ) 
dZ Z 

whence 


A r a 2 f 1 - ^ ==A T a 2 P — = ==■ 

Jo Z 2 V(l-Z 2 )(a 2 -Z 2 ) Jo V(l— Z^Ca 2 — Z 2 ) 

A r P di- 

Jo V 1-Z 2 


A7 vfb 


=NaK—NaE—N V(1 ^(a 2 P) 

The first two terms on the right are exactly canceled by t 
two terms previously obtained, so that the final result is 


i 

o 

he 


lim R.P. A 7 w 2 (z)=lim-AT -vA - ftfe 2 ^ 1 

z— >ca Z— >0 t 0 


“teWO-SX 1 -?) 


Derivation of ir 3 . — In order to cancel the effect of the 
terms Mwi+Nw 2 at infinity, the function w 3 must approach 
— R.P.A r u> 2 (z) as s increases without limit. In addition, it 
must have no singularities in the upper half of the 2 -plane, 
it must satisfy the condition of continuity at z— ± 1, it must 
be a pure real on the fixed boundaries and a pure imaginary 
on the free boundaries, and it must be zero at 2=0. It is 
readily formulated as 



That this function satisfies the first condition is readily shown 
by writing it in a slightly different form and expanding the 
radical: 



•) 


=lim— A 7 2 ■ 


N(a*-1) 


Or 


Comparison of this expression with that for lim R.PJVtt’ 2 ( 2 ) - 

Z— *C3 

shows that the difference between the two expressions 
approaches zero as z approaches infinity. That the function 
satisfies the remaining conditions is readily verified by 
inspection. 

The complex velocity function for the closed-open-closed 
tunnel with unequal pressures is, finally, 



where the modulus of the elliptic integrals is 1 ja. 

Induced velocity on the axis. — For the special case in 
which z=iy, the preceding equations for Q reduce to a 
somewhat simpler form. The procedure will be only out- 
lined herein, inasmuch as the manipulative steps are similar 
to those already described. 

Replacing z with iy in the expression for w Y and then sub- 
ct 2 

stituting y 2 =jz~a 2 reduces the first term to 


2 i 


{K'—E') 


K'-F'f— =L=\ 


The same substitutions, together with the previously described 
technique from reference 12, reduce the second term to 



The third term is found directly as 
2 iK' 


2iK' j l+y 2 

** v y i+g 


The total simplifies to the form 

Q(iy)=^\ 




( 24 ) 
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Normal velocity and normal displacement at the free 
surface. — The normal velocity on the free surface may be 
written in the following form: 



By the substitution of x 2 =a 2 —(a 2 —l)l 2 , the integrals are 
readily reduced to standard forms of incomplete elliptic 
integrals, and the equation takes the form 


■ p - gg.Hr O' (fcr)-*'®' (fcrT - 


2 iK' 



(25) 


The normal displacement, or distortion, of the free surface 
is found by integrating this expression along the free surface 
in the physical plane: 


Normal displacement at x= I-P- Q(x) d£ 

= Pl.P. Q(x)— 
Jl ' TCX 


The integral may be evaluated numerically; however, the 
third term of I. P. Q(x) is amenable to analytical treatment: 


2 iK' f* / x 2 — l 

■ a 2 , I y 


- 1 dx 2 iK' C x lx 2 — 1 


TTX 


ffiV 



1 

a 2 


; dx 




At the edge of the exit lip, where x=a, this expression reduces 
to 

NUMERICAL RESULTS 

In the following sections are described some numerical 
results that were computed by the preceding equations in 


order to show the magnitudes of the entrance and exit effects. 
With regard to sign, it will be noted that a contracloclcwise 
vortex has been considered positive, in keeping with the 
usual mathematical notation. With a positive (left-to- 
right) tunnel velocity I 7 ’, such a vortex corresponds to a 
negative lift. The corresponding induced velocity v will, 
for the tunnels under discussion, be positive (upward) ; hence, 
if e represents downwash angle, it will be negative, but the 
ratio e/Cj will be positive. 

It may be noted further that, since the complex velocity 
has been made nondimensional by dividing bj r V, the com- 
ponent v is identical with — e, in radians, and the component 
u is the fractional increase in the horizontal velocity. The 
equivalence of the two ordinate scales indicated in the plots 
of the results follows from the equation relating vortex 
strength and lift: 


.r cct 

V~ 2 

Closed-open tunnel. — In figure 17 are shown calculated 
values of the induced downwash angle along the tunnel axis 
for various positions of the lifting vortex along the axis. The 
figure shows that for £,= 1.0 and 1.5, the induced angles at 
the vortex itself (£=1.0 and 1.5, respectively) are almost 
F' 

exactly which is the value for an infinitely long open tun- 
nel; and, furthermore, the two curves are symmetrical about 
the point £=£,. In fact, within the accuracy of the plot, 
these two curves are identical with the curve for an infinitely 
long open tunnel. It may be concluded that the closed en- 
trance has no effect if the vortex is more than one tunnel 
height from the entrance. For f,=0.5, which is a more 
likely location of the wing, the induced velocity at £=£, is 
0.48r, and the curve is no longer exactly symmetrical 
about the point £=£,; however, these differences from the 
conditions for the infinitely long open tunnel are too small 
to be practically significant, so that the usual infinite-open- 
tunnel theory is still adequate for £,=0.5. For £, less Hum 
0.5, the deviations from infinite-open-tunnel theory become 
larger rapidly until, when the vortex is in the plane of the 
entrance (£i=0), the induced angle at the point £=£, is 



A similar discussion applies for the vortex in the closed 
portion of the tunnel (£i<0), although this case normally has 
no practical significance. For £,= — 1, the induced angles 
in the neighborhood of the vortex are practically identical 
with those for an infinitely long closed tunnel; however, in 
the open region (£>0), the curve is considerably different 
from that for the infinitely long closed tunnel (shown as the 
dashed curve in fig. 17). 
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Symmetrical closed-open-closed tunnel. — In figure 18 
are shown similar curves for a symmetrical closed-open- 
closed tunnel of which the length of the open section is 
1.5 times the tunnel height. All curves show a sharp reduc- 
tion in the induced angle as the closed exit is approached and 
entered; however, for £i<l, the closed exit has practically 
no effect on the induced velocities at and upstream of the 
vortex. The induced angle at the vortex decreases rapidly 
as the vortex moves downstream from about £ z =1.0 and is 
I-' 

-j- in the plane of the exit (£ t =1.5). 

Closed-open-closed tunnel with one exit lip. — Figure 19 
shows results for a tunnel similar to that Just discussed except 
that one exit lip is omitted. The two curves shown are very 
similar to the corresponding curves for the symmetrical 
condition. As was pointed out earlier, the horizontal com- 
ponent of the induced velocity on the axis is not zero for this 
unsymmetrical configuration. Values of this horizontal 
component have been plotted in figure 20 for the same two 


vortex locations as in figure 19. The values are seen to be 
very small in the forward part of the tunnel but become quite 
large in the neighborhood of the exit lip. The effect is con- 
sistent with the concept of the exit lip as a concentration of 
vortices having total strength equal and opposite to that of 
the bound vortex and serving thereby to turn the air back to 
its original direction. The fact that the two curves are prac- 
tically identical lends further support to this viewpoint. 

Comparison of the three tunnel types. — In figure 21 
are compared the induced-angle curves for £l= 0.5 and 1.0 
for the three t unn el types just discussed. It is seen that- the 
differences are slight up to about £=1.0; beyond this value the 
curves for the closed-open tunnel continue to rise, while the 
others descend rapidly. The effect of the closed,exit is some- 
what larger for the tunnel with two exit lips than for the 
tunnel with one exit lip . Although the induced angles become 
slightly negative in the downstream closed region, they even- 
tually return to zero. 





Figure 19.— Tunnel-induced angle on axis of two-dimensional closed-open-closed tunnel having one exit lip with vortex at two locations along axis. Length of lower free surface 

is 1.5 times tunnel height. 
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Figure 20. — Tunnel-induced horizontal velocity on avis of two-dimensional closed-o pen-closed tunnel having one exit lip with vortex at two locations along axis. Length of lower free 



Figure 21. — Comparison of t unne l-induced angles on axis for three types of two-dimensional tunnels. Length of open sections for the closed-open-closed tunnels is 1.5 times tunnel height 

Symmetrical closed-open-closed tunnel with unequal pres- component (or the induced angle) has an almost linear 

sures on. the two free surfaces. — By means of equation (24) variation along the axis, which corresponds to the fairly 

calculations were made of the induced vertical velocities uniform curvature of the jet that would be expected to result 

on the axis of a closed-open-closed t unn el of jet length equal from the pressure difference between the upper and lower 

to 1.5 times the tunnel height and having equal and opposite surfaces. For this same condition, the integral of the normal 

horizontal perturbation velocities ( — v. b and « s ) on the upper velocity along the free surface (equation (25)), which is the 

and lower free surfaces, respectively. The results are plot- downward displacement of the jet boundary at the exit lip, 

ted in figure 22. The curve shows that the vertical velocity was found to be 3.89tf a . 



Distance from en trance, units of tunnel height 
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Figure 22.— Ttmnel-indticed vertical velocity » on axis of symmetrical two-dimensional 
closed-open-closcd tunnel having additional velocities of -Ub and non the upper and 
lower free boundaries, respectively. Length of open section is 1.6 times tunnel height. 


For a vortex r' located at £i=0.5 in the equaJ-pressure 
case, the integral of the normal velocity along the free surface 
(equation (21)), which is proportional to the upward dis- 
placement of the jet boundary at the exit lip, was found to he 
1.20r'. Accordingly, zero displacement at the exit, cor- 
responding to the existence of a closed space above or below 
the jet (see part I), will result if the flow described in the 
preceding paragraph is superposed on the equal-pressure 

flow in such proportion that 3,89'ttj,=1.20r / ; that is, p=0.31. 

The corresponding effect on the induced angle at £ =& is 

V 

found as follows: at £=£i=0.5, pfor the equal-pressure case 


(fig. 18) is 0.48. From figure 22, —at £=0.5 for the unequal- 

pressure case is —1.44. Since 0.31 X —1.44=— 0.45, it is 
seen that, if spillage at the exit lip is prevented, the induced 
velocity in the region of the vortex is nearly eliminated. A 
similar comparison of the slopes of the curves in figures 18 
and 22 in the neighborhood of £=0.5 shows that the induced 
curvature in the region of the vortex is also nearly eliminated. 

Resume of numerical results.— The induced angle at the 
lifting vortex is essentially that for an infin ite open jet if the 
vortex is more than half the tunnel height from the entrance 


and the exit. The induced angles for case 2 (one fixed exit 
boundary) are nearly the same as for case 3 (symmetrical 
exit), so that any failure of the flow to contact the upper exit 
lip should not appreciably affect the tunnel correction. 
Finally, for case 3, if enough of the different-pressure flow is 
added to assure zero displacement of the free boundary at the 
exit (that is, if spillage at the exit is prevented, as by enclos- 
ing the space into which the spillage would normally occur), 
the induced angle at the vortex may be nearly eliminated. 


III. CIRCULAR TUNNELS 


In part III an outline is given of a general method for 
calculating the boundary effect in an open circular tunnel of 
finite jet length. The solution, involving expansions in 
Bessel functions, is somewhat similar to the solution for a 
closed circular tunnel (reference 13) but is constructed so that 
it satisfies the condition of uniformity of pressure over the 
open boundary and also the condition of continuity of 
velocity at the entrance lip. Numerical results are given for 
a lifting element on the tunnel axis. 


SYMBOLS 


£, r \ , f rectangular coordinates in units of the tunnel radius 
with origin at lifting element (see fig. 23) 

£, p, 6 cylindrical coordinates (see fig. 23) 

a, b £-coordinate of entrance and exit lips, respectively 

B variable of integration 

2 variable of integration (see reference 13) 

$o disturbance potential associated with body (or with 
vortex system) 

# tunnel-induced potential 

$ c tunnel-induced potential in closed circular tunnel 

residual potential ($— $ c ) 

J m Bessel function of the first, kind of order m 
u constant longitudinal perturbation velocity on free 
surface . 

g m U) (£) mth Fourier sine coefficient of -=r-A 




<?m®(£) mth Fourier cosine coefficient of 


h 


mn 

Usm 
Cl) 


r 


(2) 


P-1 


CD 


nth coefficient in series for g m (S) (£) 

sth zero of J m ' (not including the zero at the origin) 


mth Fourier sine coefficient of — d^o+tc) 


3£ 


pml 


mth Fourier cosine coefficient of— 


p-i 


S 

p 

L 

w 

V 

R 


tunnel-induced velocity parameter (Jy 


density of fluid 
lift of lifting element 

tunnel-induced velocity normal to the £i?-plane 
free-stream velocity 
tunnel radius 
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Figure 23. — Circular close cL-o pen-closed tunnel showing coordinate systems. 


ANALYSIS 


to find a function $(f,p, 0 ), harmonic inside the cylinder 
p=l, for which (see part I) 


d($ 0 +a>) 

dp 


=0 

fi=i 


(K®, £>b) 


c>@o+$) l _ n 
2>e p=i 


KKi) 


where the region K a , f)>& is the closed portion of the tun- 
nel and is the open portion of the tunnel. Id addi- 

tion, according to the condition of continuity at the entrance 

lip (part I) the derivative must be continuous at 

f=a. The function $ is then the velocity potential of the 
additional flow due to the tunnel boundary. 

The function $ is conveniently considered in two parts: 




INTRODUCTION 

In the analysis of the three-dimensional, circular, closed- 
open-closed tunnel, an appreciable simplification results when 
the tunnel axis lies in the plane of the horseshoe vortex. 
For off-center locations of the horseshoe vortex, or for a 
source-sink body on the axis, or for the general unsymmetrical 
disturbance, certain complications arise that are related to 
the fact that the pressure on the free boundary is then not 
equal to the pressure at ± ra in the closed parts of the tunnel. 
That is, for these cases, if the net perturbation velocity is 
zero far upstream and downstream in the closed parts of the 
tunnel, a constant longitudinal perturbation velocity n^O 
will exist on the free surface. (See parts I and II.) A 
similar complication results for a source in a completely 
closed tunnel when it is treated by the method of reference 13. 

The analysis described in the following section is applicable 
directly to the case in which the tunnel axis lies in the plane 
of the horseshoe vortex and for which the longitudinal per- 
turbation velocity on the free surface is zero. (See part I.) 
In the succee ding section are derived the additional terms 
needed for the solution of the more general problem. The 
significance of the titles of these two sections will become 
clear in the analysis. 

CYUNDRICALLY SYMMETRIC TERM OMITTED 

Boundary conditions and formal expression for $ 4 . — The 
solution is developed in cylindrical coordinates (f, p, 8 ) where 
the f-axis coincides with the tunnel axis and 8 is measured 
from the horizontal plane. The relations of these coordina tes 
to the rectangular coordinates (f, 17 , f) are indicated in fig- 
ure 23. The distance variables £, V, and p are considered 
in Tin its of the tunnel radius. 

Let (£,p, 0 ) be the disturbance velocity potential asso- 
ciated with the lifting body in unlimited space (in particular, 
the velocity potential of a horseshoe vortex). It is desired 


where is the known tunnel-induced potential for the same 
vortex system in a completely closed circular tunnel (refer- 
ence 13). The determining conditions for are then 


(1) A$ a =0 


( 2 ) 


bp 


=0 


l| 




a($oF$ c ) 

at 




(p< i) 

(£<«, f> 6 ) 

(a<K&) 


(4) 


bp 



(£=«) 


The f unc t ion may b e expressed forma fly (seereferencel3) as 
cos (26) 


where g m a) (£) and are, respectively, the rath sine and 

cosine coefficients of the Fourier series for - 5 — and g and 

c>P „-i * 

jS are variables of integration. The integrations over £ 
would, in general, have the range — <= to + 03 ; however, 
condition (2) shows that the functions g^ 03), 3 = 1, 2, are 
zero from - <= to a and from b to + 00 . The convergence 
of this function and its derivatives to the desired function 
$4 and its derivatives is discussed in the appendix of ref- 
erence 13. A modification is necessary because of the dis- 


continuity in 



that may exist at £=&. 


For this case, 


the desired convergence maybe proved for regions bounded 
away from the circle p=l, !=&. 
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The assumption of zero perturbation velocity on the sur- 
face of the jet is equivalent to the assumption that the ex- 

a@ 0 +$ c )| 


pansion of- 


af 


p=i 


in a Fourier series in 8 contains no 


term independent of 8. For this reason no m=0 term ap- 
pears in expansion (26). The next section discusses the 
somewhat special treatment that ■ is required when the 
Fourier series contains a term independent of 8. 

Evaluation of g m ^ (|). — The function 4> A given in the 
preceding equation satisfies conditions (1) and (2) regard- 
less of the precise form of the functions g m U) (|). It is now 
desired to find the functions g m {i) (£) for a<Cl<C6 such that 
will satisfy conditions (3) and (4). To this end the 
functions are represented by infinite series of the form 


“(!)=&■ 


TOO 


«) 


. it | — a 

sin rr • 


2 b—a 


-S 

n=l 


Cfl 


sin nr 


b—a 


(27) 


Since g m U) (a) thereby equals zero, condition (4) is auto- 
matically satisfied. 

Substituting this series in equation (26) gives 


m= l n=o 


t W Pmn(£,p) sin m$+h mn l2> P mn ($,p) cos md] 


where 


(28) 


and, for 7i^0, 

=£ cos J 09-6*0 

In the evaluation of these two expressions, the inner integrals 
may be found analytically and the outer integrals, which 
converge rapidly, may be found numerically. It is possible, 
however, by means of contour integration similar to that 
discussed in reference 14 to transform the infinite integrals 
into infinite series that are more convenient for the present 
purpose. The contour integration and the resulting infinite 
series are given in appendix A. 

Differentiating equation (28) with respect to l and taking 
p=l gives 


a* 


[p-l 


co co r” 

SE h mn m sin me 
m= 1 n=0 L 


ap mtt q,i) 

as 


h mn i2} cos md 


a P-nd.ir 

as _ 


(29) 


The constants h mn 05 are to be determined so thahcondition (3) 


is satisfied. For this purpose the function is 

o| |p— 1 

expanded in a Fourier series in 8: 


a($o+$ c ) 


ai 


= S [r, 

p=l m= 1 


, (I) (|) sin md+r m (i> (g) cos md] (30) 


Equating coefficients in equations (29) and (30) in order to 
satisfy condition (3) then gives 


«(S)=S h mn ™ 

n—0 


&P»,(6, 1) 
as 


0 = 1 , 2 ) 


It is assumed that the functions r m (i) (£) can be satisfactorily 
approximated by a finite number of terms of these series. 

ap.,(s,i) 


This assumption seems reasonable, inasmuch as 


as 


is bounded as n approaches infinity (see appendix B) and 
h mn U) approaches zero as n approaches infinity. Thus, 


U) 


(©• 


N 

«=0 


Cfl 


hP mn(li 1) 

ai 


The functions r m U) (|) and ^ are computed at a set 

of points {| t }, i— 0, 1, 2, ... I, where I^N. The coefficients 
h mn U) are then determined (method of least squares) so that 
the expression 


i=0 L l*=0 


aP m „(it,i)T 
ai _ 


is a minimum for all values of m and j. For each pair of 
values of m and j, this condition gives iV+1 equations for 
the iV-f 1 unknowns h m0 U) , h ml u> . . . h mN a) , as the iV+l 
partial derivatives with respect to h mn U) must be equal to 
zero. These equations are 



dPmk (Ipl) 

ai 


N 

S h « 

itf mn 

71=0 


dP mn( lt,l) dPmk (It, 1) ~| p. 

£>! at 


(k=0, 1,2, ... N) 


Remarks on the computations. — The points {|<} and the 
value of N are chosen so that the addition of more points and 
increasing the value of N will no longer appreciably affect 
the results. It is clear that the point |=6 cannot be used 
and that care must be taken not to choose too large a pro- 
portion of the points {|J in the neighborhood of f —a and 
|=6; any such attempt to describe more accurately the 
infinite values of g m (i) (!) at |=6 or of its derivatives at 
| =a and |=6 with a finite number of coefficients will cause 
a large error in the approximations to g m U) ( I) elsewhere in 
the interval a<|<6. 
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Since the functions r m (i) (£) rapidly approach zero as m ap- 
proaches infinity, the preceding equations need be solved for 
only a small number of values of m. The values of A mn c -° 
thus obtained can be used to give an approximation to the 
function (equation (28)): 

$a=S S [A m » CI> P»»(S,p) sin mff+A„,®P„»(£,p) cos md } 

m=l n=0 

(31) 

Any desired interference velocity may now be obtained by 
differentiating this series term by term and adding the re- 
sults to the corresponding interference velocity for the closed 
tunnel. 

The vertical induced velocity in the plane of symmetry 
is simply 

c& A 1 
df f=o P [0=0 

for points on the right side of the tunnel axis, or 

d—r 

for points on the left side of the tunnel axis. Inspection of 
equation (31) shows immediately that the ^-derivative of the 
second term in the bracket is zero for either case and the 
contribution of the first term is 

a) P»»&p) 
or 

— -s s (-1 rrnh mH ^p mn (i P ) 

P m=ln =0 

Furthermore, all vertical velocities on the axis itself may be 
obtained by considering only m=l, because 

lim — P mn (£, p) = 0 (m> 1) 

p-* 0 p 

The usual geometric symmetries also contribute toward 
simplifying the calculations. For example, if the horseshoe 
vortex lies in the horizontal plane of symmetry of the tunnel, 
„ un—h ®=n 

If, in addition, the vertical plane of symmetry of the tunnel 
is also the vertical plane of symmetry of the horseshoe vortex, 
all even values of m are el i mi n ated; in the corresponding 
an Asy mm etrical case (as with aileron deflection) all odd values 
of m are eliminated. 

CYLINDRICALLY SYMMETRIC TERM 

For a normal veloctiy at the tunnel wall go© that is 
independent of 8 the potential function cannot be given 


p d0 


p Oo 0=0 P m=i n*»0 


d<£ A __ J L 
df r=o p S0 


exactly in the form of the preceding section since for m=0 
the integral with respect to g will, in general, not converge 
for g in the neighborhood of zero. It is necessary to add 
additional terms to the potential so as to insure the con- 
vergence of the integrals with respect to g. Moreover, these, 
terms must be of such a nature that the potential function is 
still harmonic and gives the required normal velocity at the 
tunnel wall. 

The singularity-free potential inside the tunnel then takes 
the form 

Iwfik S '. cos 9 ® -8 <w_ 

where k— lim and where f g 0 (j 3) dp is the Cauchy 

Jfl|— P J — ® 

principal value of the integral. Both the limit and the 
integral must be assumed to exist. 

The appearance of these additional terms is not wholly 
due to the presence of the open section in the tunnel. For 
a source in a completely closed tunnel the second term does 
not vanish and would have to be used in calculating the 
tunnel-induced perturbation velocity by the method of 
reference 13. However, for a closed body or a vortex system 
plus its reflections in a completely closed tunnel, 'both of 
these additional terms vanish. 

It is easy to verify the fact that the additional terms do 
insure the convergence of the integral with respect to g. A 
straightforward differentiation then shows that $ is in fact 

harmonic and satisfies the boundary condition ^ =<7o(£)- 

For the closed-open-closed tunnel, the boundary condition 
(see part I) that the velocities far upstream and downstream 
be equal is no longer automatically satisfied by putting the 
total tangential velocity on the jet surface equal to zero. 
The determining conditions for <fu are now 


(1) 

II 

O 

(p< 1) 

(2) 

O 

II 

►-1 

II 

Q. 

io?l Q - 

lol'O 

(£<a and £>&) 

(3) 

d($o+$c) , .. 

3{ S| - _,+* 

(a^ £<&) 

(4) 

£— v-j-cc £— ► — C3 


(5) 

-0 

Op 1 

II 

.g. 


956646 — 51 S3 
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Conditions ( 1 ) and (2) arc satisfied by assuming <7o(£)=0 
for £ <a and £> 6 . Thus 

(32) 

Again it is desired to find g 0 (£) for a <£<6 so that will 
satisfy conditions (3), (4), and (5). Tbe representation of 
< 7 o (£) in tbe same form as before (equation 27) 

ffo(£)=A 0 o sin^l — -+2J«'o« smnx-r — - 

Z 0 — d n-l 0 — d 

automatically satisfies condition (5 ) . Subs tituting this series 
in equation (32) gives 


There exists a unique value of u for which the coefficients 
h 0n will satisfy this equation and it can be found as follows: 
Let 

h 0 „=h 0n '+uh 0n " 


where 


co 


sv 

71 = 9 


ajVfe i) 

a? 




a? 


p-i 


and 


Then if 


S Kn" 
7»=0 


aPon(^,i) 

a? 


i 


h m '+uh 00 " . ^ (Aon' +uh on ")[l — (— 1)"] __ Q 

X ’nti 7 lx 

2(6— a) 6— a 


— S honPan (£j p). 


7i=0 


where 


P»(l, ri-i X’ ‘ s ” I fci cos «« “0 «“ 

P..«, «)-i J 0 “ [A|^j£ S in.x cos <tU-0de- 


2 r 6 


smnx 


P~ 


^ d/3 dg 


The resulting infinite series for the P 0 n(£,p) obtained by 
means of contour integrations are given in appendix A. 
Condition (3) then becomes 


co 


2}A 

71=0 


Ore 


dJWM) _ 

d($o+$c) 

a? 

d£ 


+U 




But (see appendix A) 


Iim 

£-»-{- a> 



lim 

£— > — co 



and so 


r a<fu 

hm lim ~ 


so that condition (4) becomes lim ^^= 0 . Thus it is nec- 

|-*+® 0 £ 

essary that 


__Aoo__j_y^ jUl ( 1) ”] q 

7 T 72=1 7l7T 

2(6— a) 6— a 


u 


2Aoo'+S ~~ [1 — ( — 1)“] 

71=1 n 

2 Aoo" + S ““ [1 — ( — 1)”] 
Ti-i n 


The coefficients A 0 / and k an " are found by solution of 
sets of simultaneous linear equations, as described in the 
previous section. 

The function = Y^.hcJ'PnJL p) is the perturbation 
potential which, when added to that of a uniform flow, gives 
the potential of the disturbance-free expanding tunnel 
described in part I and indicated in figure 2 . The cor- 


responding perturbation velocities 
opposite values at « and — » . 




have equal 


and 


EXAMPLE 


As a somewhat simplified illustration, the problem of a 
semi -infini te unit doublet distribution (degenerate horseshoe 
vortex) along the tunnel axis was considered. The tunnel 
was assumed to have an open jet 3 tunnel radii in length. 
The tunnel interference was calculated for four different 
positions of the upstream end of the doublet distribution, 
these positions being 0.1, 0.4, 0.7, and 1.0 radius downstream 
from the entrance. If the upstream end is taken as the 
origin of the coordinate system, then (see appendix C) 


r '“’® — (MlT’-ir ZJ?$i ) C08 « f 

( 33 ) 

where K 0 and K x are the Bessel functions as defined in ref- 
erence 15. 

The points {£<} were taken as a, a+0.3, a+0.6, . . . 
<z-{- 2.7, that is, a set of tenpoints, at 0.3 intervals, starting with 
the entrance lip of the tunnel. The coefficients Ai„ (1) were 
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Distance, from lifting element , units of tunnel radius 


Figure 24. — Turmel-Induced velocity parameter along tunnel axis for several positions of the lifting element in a close d-open-closed circular tunnel. 


found by the method of least squares for A T =0, l, , 5 
and also so as to satisfy the equations 

71 = 0 0 ? 

at all ten points. Plots of the resulting functions gi m (Jb for 
the different values of N indicated that convergence was es- 
sentially complete for N between 3 and 5. This simplifica- 
tion results in appreciable saving in the amount of compu- 
tation. Not only is it necessary to solve a smaller set of 
simultaneous equations but also P ln and bP lre /c)t need be 
found for fewer values of n. 

The computation was fairly straightforward. In the de- 
termination of K 0 and Ki were obtained from the 

tables of reference 15, and Ji and J-l were obtained from 
the tables and from the relations between the Bessel func- 
tions and their derivatives (references 15 and 16). Weddle’s 
formula (reference 17) is convenient for performing the inte- 
grations. In the case of Pm, the values of y 3l appearing in 


the formula for <2i» w (p) were found from the formula in ap- 
pendix III of reference 15, and Ji and J/, as just noted. 

In the evaluation of Iim — P ln it is noted that the value of 
P-0 P 

um 1^=1. 

The results of these computations together with those for 
the completely open and completely closed tunnels and those 
given in reference 3 are shown in figures 24, 25, and 26. In 
figure 24, the vertical tunnel-induced velocity along the axis 
for the four different positions of the lifting element together 
with the results for the open and closed tunnels are plotted 
against distance from the lifting element. The same results 
are plotted against the longitudinal distance from the 
entrance lip in figure 25. Figure 26 shows the results of 
reference 3 compared with the results of the present report 
for the same case — that of the lifting element 1 radius 
downstream from the entrance lip. 

The tunnel-induced velocity in the upstream regions and 
in the neighborhood of the lifting element, although only 
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Fioure 25.— Tunnel-Induced velocity parameter along tunnel axis for several positions of the lifting element in a closed-open-closed circular tunnol. (Same curves as in fig. 21 but 

plotted against distance from entrance.) 


slightly less than that for an open tunnel for the lifting ele- 
ment 1 radius downstream, falls off more and more rapidly 
as the lifting element is moved towards the entrance lip. 
The maximum induced velocity is attained about 1 radius up- 
stream from the exit and is never more than 78 percent of 
the maximum value for a completely open tunnel. After 
the maximum the values fall rapidly and approach the 
values for a closed tunnel in the downstream regions. The 
results of reference 3 (see fig. 26) are consistently below the 
present results especially in the region behind the lifting 
element. 

An extrapolation from the present results indicates that 
the induced upflow at the lifting element, for the lifting ele- 
ment in the plane of the entrance lip, is approximately zero, 
or the average between the completely open and com- 
pletely closed cases. The same result (that the effect in the 


plane of the entrance lip is the average of the effects for the 
completely open and the completely closed tunnels) was also 
obtained for the two-dimensional tunnel (fig. 18). 

CONCLUSION 

For an open wind tunnel, the corrections corresponding to 
an infinitely long open jet will usually be adequately accurate 
if the region of interest (where the lift is located and where 
the boundary-induced flow is being considered) is at least 
half the jet height from the jet entrance and exit sections. 
As the distance of the lifting element from the entrance is 
decreased below this limit, the boundaiy-induccd flow 
decreases rapidly and, when the lifting element is in the 
entrance plane, the induced angle at the lifting element is 
about the average of that for an open tunnel and that for a 
closed tunnel. 
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Figcbe 26 —Comparison of the results o£ the present report with those of reference 3 for a= — 1 
together with those for open end closed circular tunnels. 


In the theoretical studies of these flows, the usual boundary 
conditions of pressure uniformity on the free surface and of 
zero normal velocity on the closed surface must be supple- 
mented with the conditions that the velocity be continuous 
at the entrance lip and that the velocities far downstream and 
far upstream in the closed sections be equal. For the two- 
dimensional open tunnel, a convenient general mathematical 
approach is to transform the infinite strip (representing the 
tunnel) to the upper half-plane by the logarithmic trans- 
formation and then to develop the desired complex velocity 
in this transformed plane. For the circular open tunnel the 
solution may be effected by expressing the potential by a 
finite series of Bessel functions, satisfying the boundary condi- 
tion on the free surface at a finite number of points, and solv- 
ing for the coefficients by simultaneous linear equations. 

For noncircular open wind tunnels, solutions in terms of 
available functions will be very inconvenient. For such 
cases, the trends indicated by the present results may suffice, 
when applied to the presumably known corrections for the 
in fini tely long open and closed configurations, to provide 
adequate corrections. Solutions for the general three- 
dimensional configuration may also be possible by electrical- 
analogy methods, in which either the perturbation velocity 
potential or the acceleration potential is analogous to the 
electrical potential in an electrolyte solution. Such analogies 
may be characterized, however, by considerable technical 
difficulty. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 20, 1948. 


APPENDIX A 

EVALUATION OF P m ,(f, p) 


EVALUATION FOR m^O 

Contour integration. — For n^O, 

J mij'PQ) 




3 (I, 

sm mrj cos g(/3 — 0djS 


The inner integral may be found directly: 


(Al) 


6 jS — d 

sinner cos g(|8 — |)d/3 


b—a 


i, 


UK 

b—a 


„2_( 

nx V 

2 k 

h-a) 


[(—!)* cos 2(5 — 0-cos g(£— a)] 


(A2) 


The problem of evaluating P mn (i,p) thus reduces to that of 
evaluating integrals of the form 


C^ Jmjipq) cos kqdq 

Jo i( 


iqtJm'iitiW-h 2 ) 
Consider the integral in the complex 2 -plane 

d m (ipz)e m 


1 T 


2%iJ izj m '(iz)(z 2 -h*) 


dz 


around the contour indicated in figure 27. Its value is the 
sum of the residues of the integrand at its poles inside the. 
contour. These poles are the values of z for which J m ' (is) = 0. 
The zeros of JJ will be designated y„ a ; they are real and 
may be obtained from the formula in appendix III of refer- 
ence 15. The poles of the integrand then occur at iz=y sm , 
that is, at z= —iy sm - Since only the poles within the contour 
are desired, only the negative zeros of d m ' are considered. 

The residue of — — — at s= —iy, m is 


lim 


Z + iVsr, 


lim 


z+iy , 


iy*m d vi (iz) z-*—iy, m d m ( is ) dm'tysm) 

since J m '(y sm )= 0; by the definition of the derivative this 
expression reduces to 


n) 



The residue of the integrand at s= —iy, m is thus 

■ Jm(py>d)e kv ‘ m 
y,m(y sm 2 +ld)J m "(y sm ) 

But the Bessel functions satisfy the relation 

•v'+j * r -'+( 1 “$) J '-“0 

so that at x=y sm , where J m '= 0, 


jm —y, m 
y* m 2 




dm(y sm) 


d vi ( p y e k ' J ‘vy s 


d 71 (i/jjb) 

whence, finally, 

1 yfi d m (ip z)e ila dz _ 

27 ri J izd m , (iz)(z 2 —h 2 )~ 1 ' V (y sm 2 +h*) (m 2 -y im *) d m (y, m ) 

where the number of terms in the summation depends on 
the radius of the outer semicircle. 

For k 0, if the radius of the outer semicircle is allowed to 
approach infinity in discrete steps so as to avoid the poles of 
the integrand, the integral over the semicircle approaches 
zero. The limiting values of the integrals along the two 
inner semicircles, as -their radii approach zero, are readily 
determined by the usual process as 


and 


• d m (—iph)e ikh 
U 2 d m '(-ih) 

■ d m (iph)e ikh 
4:h 2 J m '(ih) 


540 



LINEAR THEORY OF BOUNDARY EFFECTS IN OPEN WIND TUNNELS. WITH FINITE JET LENGTHS 


541 


These two terms may be combined after noting that reversing the sign of the argument in J m and JJ merely reverses the 
sign of their ratio, to 


J rr. ( 1 P IQ 

2h 2 JJ(ih) 


sin kh 


Equating the total integral along the infinite contour to the sum of the residues thus gives 

1 f“ J m (ipq)e m dq JJiph) ■ M ^ J m (j>y,J)e.~ t *™y,m 
2 Hj- a iqJJ(iq)(.q>-h*)~2h>JJ(ih) ^ U (j/.J+h^rrr-y^JJy^) 

where the y sm terms are now defined as the positive zeros of JJ instead of the negative zeros (if JJ(x)=Q, so also does 
JJ(—x)). By equating the imaginary part of the left-hand term to the right-hand side, which is a pure imaginary, 
there results, finally 

J n (ipq) coskg dq /iriJ m (iph) . ,, ^ JJpysJe-^y^ 

Jo iqJJ(iq)(q*-h?) 2h*JJ(ih) sm (y,J+h*)(;m*-y,J)J a (.y m ) 

Expressions for P mn , 0, m?^0. — In the preceding development it was assumed that k^O, which was acceptable with 
regard to equation (A2) in view of the fact that the cosine is an even function of the variable. This essentially non- 
negative value of k must be retained, however, in the final expressions for P mre (£, p): 


P =_ i 

x mn cy 
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. fix T , / . Tlx \ 

1 J “ (* b=s) 
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b—a b 
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Now for a^^b. 


(A3) 


For Ka, 


For£>6, 


(-!)• smfcfe=(-l)‘ sin 
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b — a 


||-a| 

= — SlUTlx ^ L 
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a — ? 
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\I=A 

b—a 
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§ — a 


=sm nx 


b—a 


It-a| 

=sm nx 1 

b—a 
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The first term on the right-hand side of equation (A3) is thus equal to zero for £<Ca or £>6. The desired expressions^ for 
Pmnii, p),n^0, are, therefore, for a^^b, 


Pmn(%) P ) — 


T / . mr \ 

Tm \ %p b-a) 

. mr , mr \ 

'T=i J -X % T=i) 


o 


sinwTrl— 

0 Ur 0 Or s~0 


and, for |^a or |^J, 


where 


p)=- 7 r 5 --S 

W 5=0 

<2»» (S) =P , ^(py sm )y im ( ft j^ 0) 


J(m 2 — y,^ 2 ) 

The calculations for n= 0, which follow similar lines, are not given here. The final formulas are for a sji £<(6, 


Pmo(£t p) — 


‘0 P 2(6-a)) 


for 


and for £>6, 


where 



f i/it \ 

2(6— a) 1 

\2{b—a)J 
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t 2{b-a) e 


si ” f H+S C-»w[r-.— 


<2mO (S) (p)==- 


Jm(py sm)y i 


| (® 2 ysm^Jmiys 


EVALUATION FOR m=0 


The evaluation of P mn (£,p) for m=0 proceeds essentially as before with the difference that the contour of integration 
must avoid the origin. For n= 0, the final formulas are as follows: for aSs £<(6, 


for £<«, 


f— <* 

j 1/0 [" P 2(6— a) J 

7T 

. r/ ,• v 

2(6— a) 

2(6 -a) ^ L 2(6 - a) J 


sm 


2 & 


=f+§ «*•" «-«-»•«] 


FV s j+-y-+a «.» 

2(6— a) 
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and for £>5, 




2(5— a) 

For nr* 0, the corresponding formulas are as follows: for ag £<(5, 


P Ottfep) = 


for 


and for £>6, 


j ( i 7117 ^ 

1 [tt-a)-(-l)*(&-e] + — L-- ■ ' ^ — ^sinKrl-—— t~~S 

/ / * W 7 T \ 0 — & 0 — ff s =0 


TIT 
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. n 7T T,/ - 71 IT \ 


- ?i t o — a »=o 

6— a 

Ponfep)= — [(*-«)-(- 6 ) 1 - S &.« [e- (— i)» e 

7ix o — a c=o 

6 — a 


APPENDIX B 

PROOF THAT is BOUNDED AS w-*a>,a<|<b 


Differentiating the formula from appendix A for P mjI (!,p) and putting p = l gives 


eos nic J?fL_ 2 ^ 

o— a o— a S 




w ( 4 wj) 


^jf.o , +( 6 _ a ) , J (m ,! y*»*) 


[g-C£-«y„^(_ !)* g-C6-£)y„„] 


The second term of the right-hand member is of order— for large n so that 

n 


/ . 71X 

y b—a 


Urn ^f^U ixn - 

n— »<n Oq n— *•«= • t r( • TIT ^ 

V 6-o/ 


COS71TT 


£— a 
b — a 


The cosine factor of this expression merely oscillates between 1 and —1. For the remaining factor, it is noted from the 
asymptotic expressions for the Bessel functions (reference 15, pp. 59-61) that is essentially of the form i m —== as 

f-s-m , from which it can he readily shown that the fraction approaches unity as t—> ®> . 

i Jm V-*) 
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APPENDIX C 

DERIVATION OP EQUATION (33) 


Equation (33) was derived for use in calculating some of 
the results given in reference 13, but it was not explicitly- 
stated and discussed in that paper. Because certain steps in 
its development are not obvious, the present outline of its 
derivation is given. Familiarity with reference 13 will be 
assumed. 

Certain difficulties arise in the treatment of the doublet 
line directly; so the result is found by considering a horseshoe 
vortex of finite span and letting the span approach zero. 

Equation (6) of reference 13 gives the formula for 


bp 


p=i 


(where $2 is defined in reference 13) corresponding to a horse- 
shoe vortex of strength P and span <r having one trailing 
vortex along the tunnel axis and the other to the right of 
the axis. The procedure for the doublet consists of letting 
the yaw angle rp be zero, expanding the radicals in ascending 
powers of <r, and proceeding to the next step in the analysis, 
where <r will eventually be made to approach zero. In the 
expansions, powers of er higher than the first may be neglected 
except where <r occurs in the product £<r, since £ takes on 
infinite values; furthermore, since for the doublet the field 
should be symmetrical about the vertical plane of symmetry 

unsymmetrieal factors, as & cos 0, may be imme- 


diately eliminated. 


The formula for ^ 

dp 


is thus 


P = 1 


dS> 2 

bp 


P = 1 


==— E^lim (sin d 

,_»g ( 


( 




1 


Vl + ?V Vl + ! ! 


■> 


£ sin 0 / 1 \ 

<r(£ 2 +sm 2 8 )\ VT+£V ) 

| sin 0 p 1 cos 2 0 1) 

s 2 +sm 2 0Lvr+^ a + a s/2 j 


(Cl) 


of equation (8) of reference 13 gives 

S/-' (v^p-vfTp) 0085 ®- 8 

which, if integrated by parts, reduces to 

2 J— [(1 +^ 2 ) 3/2— (1 + / 3 2 <r 2 ) 3/2 ] SU1 d P 

The contribution of the first term in the brackets is 

If” singQ3-|) dp „ . 

gj_„ (1 + /3 2 ) 3 ' 2 ~ 2 ^( 2 ) singj 

(See reference 15, p. 52.) The fact that the contribution of 
the second term in the brackets is zero follows immediately, 
upon performing the change of variable p—P<r, from the 
Riemann-Lebesgue lemma (reference 11, p. 172). 

The third term in the braces of equation (Cl) is converted 
as follows: 

£ sin 0 I - 1 cos 2 0 ~] £ sin 0 

£ 2 +sm 2 0 L vr+y 2 a + a 3/2 J (IW ! 

so that again the first and only Fourier coefficient is given 
directly. Inserting it into the inner integral of equation (8) 
of reference 13 and integrating by parts gives 

= — 2q[K 0 (q) sin g£ 

The second term in the braces of equation (Cl) is not a 
one-term Fourier series. Its nth Fourier coefficient is 
given by a constant times 


According to the procedure of reference 13, it is necessary to 
make a Fourier analysis of the three terms in the braces and 
then to insert the Fourier coefficients in equation (8) of 
reference 13. 

In the first term in the braces, the expression within the 
parentheses is the first and only Fourier coefficient. Chang- 
ing £ to j3 and inserting the expression into the inner integral 


C 2 * £ sin 0 sin nd / 1 \ , 

Jo o-(£ 2 +sin 2 0) \ yl + £V/ 

Inserting this expression in the inner integral of equation (8) 
of reference 13 and reversing the order of integration gives 
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After substitution of p=P<r, the limit of the inner integral 
becomes 


lim 

TWO 


\ J". p+#3p*( i_ vrh) cos 2 (H ) iv 


Integration by parts and elimination of terms in <r reduces 
this expression to 


lim— 

<r-H) 2 



2p*+l 1 • 

p 2(l_L^2)3/2j S1Ii 2 




which is zero, by the Riemann-Lebesgue lemma. 
Einally, then, for the unit doublet 


Jifei ) [2Kt(2)+ ™ 1 ™ 2Si2 

The potential $ 0 of a unit doublet line along the axis is 




““r+^W^+^+r 2 

whence 


/ £ i 1 N \_sinfl/__i_ I A 

\VP+n , +r , " r ) P \Vl 2 +P 2 ) 


sin 0 

"a? t=rO; 2 +i) 3/2 


The flow of the usual reflection vortices for wings of finite 
span reduces, as the span becomes arbitrarily small, to a 
uniform upflow in the finite section of the tunnel and there- 
fore contributes nothing to the longitudinal velocity. 

The coefficient of sin 6 in — is thus seen to 

p= i 

be the expression given in equation (33) . 
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